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II.  ABSTRACT 

The  important  wave  parameters  to  be  considered  for  tsunami  model  studies  are  wave 
height  and  period,  and  wave-front  orientation.  The  first  two  of  these  parameters  can 
be  determined  by  marigraphic  measurements  or  by  visual  observations;  however,  wave¬ 
front  orientation  has  never  been  accurately  observed  at  the  problem  site  (Crescent  City, 
Calif.).  A  digital  computer  program  was  written  to  plot  wave  rays  from  three  recent 
epicentral  locations  to  Crescent  City  to  obtain  approximate  tsunami-front  orientations. 
The  refraction  diagrams  were  checked  by  comparing  the  computed  and  actual  arrival  times 
of  the  wave  fronts.  The  actual  arrival  times  were  obtained  from  recording  tide  sta¬ 
tions  at  Crescent  City  and  iftlo,  Hawaii.  Initial  wave-front  orientations  near  the 
earthquake  epicenter  were  either  assumed  or  taken  from  the  literature  The  frequency 
distribution  of  wave  heights  at  Crescent  City  is  a  necessary  consideration  in  the  over¬ 
all  study.  The  period  of  record  for  tsunamis  at  Crescent  City  is  too  short  to  deter¬ 
mine  the  frequency  distribution  of  the  wave  heights;  thus,  an  attempt.-was  made  to  cor¬ 
relate  the  tsunami  records  from  particular  stations  in  the  San  Francisco  and  Hilo  Har¬ 
bor  areas  so  that  the  Crescent  City  record  could  ~he  artifiM-ffTTv  lengthened..  Correla¬ 
tion  between  the  Hilo  Harbor  and  C r e s c enJ^City-datffTndic ate s  that  the  1964  tsunami  at 
Crescent  City  was  abnormall5H§(?vere .  Effects  of  distance  on  wave-height  attenuation 
were  invest igated^and  compared  with  theoretical  predictions.  On  the  basis  of  these 
considerations^  a  tentative  frequency  relation  was  derived  and  a  risk-duration  relation 
was  prepared  from  the  frequency  relation.  The  selected  values  of  the  test-wave  dimen¬ 
sions  (height  and  period)  and  orientation  determine,  to  a  considerable  extent,  the 
model  configuration  necessary  to  ensure  results  sufficiently  accurate  for  purposes  of 
the  model  study.  Since  experience  in  the  design  and  operation  of  long-period  wave 
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ABSTRACT  (Continued) 

models  is  scant,  it  was  deemed  necessary  to  conduct  preliminary  tests  utilising  a 
pneumatic  wave  generator  in  a  two-dimensional  flume.  The  results  of  these  tests 
were  compered  with  theoretical  predictions  from  an  idealized  matheiAfcical  model 
and  on  the  basis  of  these  tests,  a  preliminary  model  design  was  proposed.  Further 
two-dimensional  testing  is  deemed  necessary  before  a  satisfactory  design  of  \,ue 
three-dimensional  model  can  be  accomplished.  _ 
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The  theoretical  study  described  herein  was  requested  by  the  District 
Engineer,  U.  S.  Army  Engineer  District,  San  Francisco,  in  a  letter  to  the 
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indorsement  dated  7  November  1967*  The  investigation  was  conducted  during 
the  period  December  1967  to  August  3.968. 
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a  Radius  of  zero  contour  of  the  theoretical  island  in  fig.  Cl; 
constant  used  in  an  integration  formula 
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A  Constant  used  in  equation  Dl? 


Coefficients  used  in  the  surface  fitting  procedure  of  the 
computer  program 


2  J 


b  Radius  of  the  largest  contour  of  the  theoretical  island  in 
fig.  Cl;  constant  used  in  an  integration  formula 
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Physical  distance  between  wave  rays  at  a  point  of  unknown  wave 
height,  h^ 

Physical  distance  between  wave  rays  at  a  point  of  known  wave 
height,  hQ 

Constant  used  in  equation  Dl4 
Effective  bandwidth  of  band-pass  filter 
Tsunami  celerity  in  the  open  ocean;  C  - 


Constunt  found  in  equation  B9;  C^ 


Cg  Constant  found  in  equation  C15 

D  Time  in  years  in  equation  2 6 

e  Probable  error  inherent  in  the  estimate 
8 

f^  Fundamental  frequency  of  a  resonating  body;  equation  1 

h  height  of  initial  wave  from  generator;  elevation  above  mean 
water  level 

h^  Unknown  wave  height 

hc  Open-ocean  wave  height  near  Crescent  City  shelf 

^  Open-ocean  wave  height  near  Hilo  shelf 

hQ  Known  wave  height 

h^.  Wave  height  after  t  hours  of  travel 


ix 


Wave  height  after  one  hour  of  travel 

Depth  of  water  through  which  wave  ray  is  passing;  assumed  water 
depth  on  a  Mercator  map;  mean  water  depth  at  any  location  in  a 
channel 

Depth  value  assigned  to  an  intersection  within  a  grid  of  depth 
values 
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K 
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*1 

l 

L 

m 

n 
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Hx,y) 


Initial  water  level  in  flume  before  release;  water  depth  at 
wave  generator;  water  depth;  water  depth  over  a  shelf  at  a 
distance  L  from  the  shore  (equation  Dl) 

Water  depth  at  a  random  point  within  a  grid  of  depth  values 

Appointed  depth  defined  in  fig.  Dl;  maximum  depth  of  ocean 
surrounding  theoretical  island 

Bessel  function  of  order  zero  of  the  first  kind 

Bessel  function  of  order  one  of  the  first  kind 

Constant  between  0.06  and  0.10  depending  on  the  air.  space  and 
the  aperture  openings  of  the  pneumatic  wave  generator 


C  2-2 

Constant;  K  =  ^  =  ^1  “  >  radius  of  curvature  of  a 

H  A2 

curve  in  rectangular  coordinates;  constant;  K  = 

Bessel  function  of  order  zero  of  the  second  kind 
Bessel  function  of  order  one  of  the  second  kind 
Length  of  sloping  channel 

Water-surface  length;  length  of  submerged  shelf  off  Crescent 
City 


Curve  fitting  coefficient 

Curve  fitting  coefficient;  ratio  of 
freedom 


tx"  311(1  5  degree  of 


Average  frequency  of  occurrence  of  tsunamis  of  equal  or  greater 
height  h  ;  factor  used  in  the  solution  of  equation  D30 


N  = 


jfo  _ Utt 
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Arbitrary  point  or.  a  curve  in  x  ,  y  coordinate  system 
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Successive  computation  points  along  a  wave  ray 
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q  Risk  associated  with  an  equal  or  greater  event  occurring  during 
a  given  time  period  D  years 

r  Polar  coordinate  measure  of  distance  r  =  ■— ;  radius  of 
hypothetical  island  between  zero  and  maximum  contour 
.  .  dr 
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r  r  =  — g 

ae 

R  Radius  of  the  earth,  assuming  spherical  shape 
s  Distance  along  a  wave  ray 
S  Vertical  model  scale 
SL  Sample  length 

S,  Ray  path  over  a  given  Mercator  map 


t  Time;  tsunami  travel  time  t 
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~  Js-  V 


T  Period  of  forcing  function;  wave  period 

T  Fundamental  period  of  bay  oscillation;  tsunami  travel  time  from 
source  to  the  Crescent  City  area 
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cr 

Th  Tsunami  travel  time  from  source  to  Hilo  area 
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x  Measurement  of  length;  part  of  total  length  of  water  surface  L 
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z  =  sinh  y 
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Dimensionless  depth  parameter 

Coordinates  of  the  random  point  measured  from 
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upon  opening  air  valves 
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equation  D12 
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British  units  of 
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O.092903 

0.3048 


converted  to  metric 

To  Obtain 

centimeters 

meters 

kilometers 

square  meters 

meters  per  second 
per  second 
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SUMMARY 


The  important  wave  parameter/3  to  be  considered  for  tsunami  model 
studies  are  wave  height  and  period,  wave-front  orientation,  and  the  fre¬ 
quency  of  occurrence  of  waves  of  different  heights.  The  first  two  of  these 
parameters  can  be  determined  by  marigraphic  measurements  or  by  visual  ob¬ 
servations;  however,  wave-front  orientation,  which  is  an  important  variable 
in  the  construction  and  operation  of  a  model,  has  never  been  accurately 
observed  at  the  problem  site  (Crescent  City,  Calif.).  A  digital  computer 
program  was  written  to  plot  wave  rays  from  three  recent  epicentral  loca¬ 
tions  to  Crescent  City  to  obtain  approximate  tsunami-front  orientations 
and  to  compute  refraction  coefficients  in  -the  Crescent  City  area.  The 
refraction  diagrams  were  checked  by  comparing  the  computed  and  actual  ar¬ 
rival  times  of  the  wave  fronts.  The  actual  arrival  times  were  obtained 
from  recording  tide  stations  at  Crescent  City  and  Hilo,  Hawaii.  Initial 
wave-front  orientations  near  the  earthquake  epicenter  were  either  assumed 
or  taken  from  the  literature. 

The  frequency  distribution  of  wave  heights  at  Crescent  City  is  a 
necessary  consideration  in  the  overall  study.  The  period  of  record  for 
tsunamis  at  Crescent  City  is  too  short  to  determine  the  frequency  distribu¬ 
tion  of  the  wave  heights;  thus,  an  attempt  was  made  to  correlate  the  tsu¬ 
nami  records  from  particular  stations  in  the  San  Francisco  and  Hilo  Harbor 
areas  so  that  the  Crescent  City  record  could  be  artificially  lengthened. 
Correlation  between  the  Hilo  Harbor  and  Crescent  City  data  indicates  that 
the  196U  tsunami  at  Crescent  City  was  abnormally  severe.  Effects  of  dis-  ' 
tance  on  wave -height  attenuation  were  investigated  and  compared  with  theo¬ 
retical  predictions.  On  the  basis  of  these  considerations,  a  tentative 
frequency  relation  was  derived  and  a  risk-duration  relation  was  prepared 
from  the  frequency  relation. 

The  selected  values  of  the  test-wave  dimensions  (height  and  period) 
and  orientation  determine,  to  a  considerable  extent,  the  model  configura¬ 
tion  necessary  to  ensure  results  sufficiently  accurate  for  purposes  of  the 
model  study.  Since  experience  in  the  design  and  operation  of  long-period 
wave  models  is  scant,  it  was  deemed  necessary  to  conduct  preliminary  tests 
utilizing  a  pneumatic  wave  generator  in  a  two-dimensional  flume.  The  re¬ 
sults  of  these  tests  were  compared  with  theoretical  predictions  from  an 
idealized  mathematical  model,  and  on  the  basis  of  these  tests,  a  prelimi¬ 
nary  model  design  was  proposed.  Further  two-dimensional  testing  is  deemed 
necessary  before  a  satisfactory  design  of  the  three-dimensional  model  can 
be  accomplished. 


xv 


THEORETICS  IN  DESIGN  OF  THE  PROPOSED 


CRESCENT  CITY  HARBOR  TSUNAMI  MODEL 


PART  I:  INTRODUCTION 

1.  Crescent  City,  Calif. ,  with  a  resident  population  of  approxi¬ 
mately  3200,  is  located  on  the  northern  California  coast  approximately 

20  miles*  from  the  Calif ornia-Oregon  border.  The  combination  of  nearshore 
undersea  topography,  resonant  characteristics  of  the  surrounding  nearshore 
area,  and  exposed  position  on  the  coast  makes  Crescent  City  particularly 
susceptible  to  tsunamis  originating  anywhere  in  the  Pacific  seismic  belt. 

2.  The  primary  purposes  of  this  investigation  were  to  conduct  a 
theoretical  study  of  tsunamis  with  particular  emphasis  on  their  occurrence 
at  Crescent  City,  and  to  provide  input  data  for  a  hydraulic  model  which 

is  to  be  used  to  determine  the  characteristics  of  proposed  barrier  plans 
to  protect  the  city  and  harbor  of  Crescent  City  from  future  tsunamis. 
Because  of  the  magnitude  of  this  many-faceted  problem,  it  was  decided  to 
divide  the  study  into  four  major  sections,  namely,  the  wave  input  for  the 
proposed  model,  the  tsunami  refraction  problem,  the  frequency  of  occurrence 
of  tsunamis  of  various  heights  at  Crescent  City,  and  the  flume  tests  con¬ 
ducted  to  determine  details  of  the  proposed  model. 

3.  In  the  Crescent  City  Harbor  tsunami  model,  as  in  any  model  of 
the  same  scope  and  nature,  one  of  the  primary  elements  to  consider  is  that 
of  wave  input.  In  a  general  sense  this  involves  decisions  as  to  the  mag¬ 
nitude,  form,  and  period  of  the  critical  waves  in  a  specific  tsunami  wave 
train.  Equally  important  is  the  orientation  of  tho  wave  fronts  approaching 
the  harbor  area,  which  determines  the  alignment  of  the  wave  generator  dur¬ 
ing  model  testing.  Unfortunately,  owing  to  the  early  destruction  of  the 
harbor  tide  gages  during  the  1964  tsunami,  the  magnitude  and  form  of  the 
waves  for  this  occurrence  could  not  be  accurately  determined.  The  disas¬ 
trous  effects  of  this  particular  tsunami  appear  to  have  been  associated 


*  A  table  of  factors  for  converting  British  units  of  measurement  to  metric 
units  is  given  on  page  xiii. 
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with  extensive  inundation  by  a  rising  horizontal  water  surface  instead  of 

1 

a  borelike  configuration  sometimes  associated  with  tsunamis.  The  limit  of 
the  inundation  is  well  defined,  and  accordingly,  verification  of  the  model 
processes  will  rely  solely  on  the  limits  of  this  inundation. 

4.  To  determine  the  wave  directions  at  the  positions  of  the  model 

wave  generator,  a  new  tsunami  refraction  diagram  technique  was  developed  by 

the  authors.  Refraction  diagrams  were  prepared  for  three  recent  tsunamis 

(Alaska,  1964;  Chile,  i960;  and  Kamchatka,  1952),  and  wave-front  patterns 

charted  from  each  of  the  sources  were  determined  for  the  nearshore  areas 

at  Crescent  City  and  Hilo,  Hawaii.  Van  Dorn's  10-mln  wave- front  configura- 
2 

tion  was  used  as  the  starting  position  for  ray  trajectories  for  the  Alas¬ 
kan  tsunami .  Starting  lines  for  the  other  two  tsunamis  were  assumed  by  the 
investigators . 

5*  In  the  selection  of  test  waves  for  the  Crescent  City  Harbor 
tsunami  model,  the  question  arises  as  to  whether  the  effectiveness  of  pro¬ 
tective  structures  should  be  determined  for  tsunami  intensities  greater 
than  the  one  that  occurred  in  1964.  In  this  connection,  a  study  was  made 
to  ascertain  the  probability  of  the  occurrence  of  tsunamis  at  Crescent  City 
with  more  critical  dimensions  than  those  of  the  1964  tsunami.  The  fre¬ 
quency  distribution  of  tsunamis  at  Hilo  Harbor  was  determined  from  a  con¬ 
siderable  number  of  occurrences,  and  it  was  assumed  that  the  Hilo  fre¬ 
quency  distribution  is  similar  to  that  for  Crescent  City.  The  correlation 
of  the  observed  tsunami  heights  for  these  two  localities  was  examined  in 
light  of  the  tsunami  sources. 

6.  After  the  orientation,  height,  and  period  of  the  test  waves  have 
been  selected,  the  linear  scales  of  the  model  and  the  extent  of  the  sea 
area  lying  between  the  wave  generator  and  the  harbor,  which  must  be  in¬ 
cluded  within  the  model  limits,  can  be  determined.  The  extent  of  the  re¬ 
quired  sea  area  is  related  to  the  number  of  wave  cycles  to  be  simulated  In 
the  operation  of  the  model.  The  greater  the  sea  area,  the  larger  the  num¬ 
ber  of  cycles  that  can  be  reproduced  in  the  model,  and  the  more  accurate 
will  be  the  observable  wave  phenomena.  Cost  considerations,  however,  Im¬ 
pose  restrictions  on  the  allowable  space;  yet  it  is  essential  that  the 
model  sea  area  be  of  such  size  that  a  minimum  of  one  initial  wave  cycle 


is  completed.  Therefore ,  the  necessary  condition  is  that  the  fundamental 
period  of  the  resonant  oscillation  of  the  model  sea,  between  the  coast 
and  the  generator,  correspond  to  the  largest  period  of  the  prototype  waves 
a elected  for  testing.  Tests  in  a  two-dimensional  wave  flume  became  nec¬ 
essary  to  develop  the  basic  ideas  concerning  this  aspect  in  the  design  of 
the  proposed  three-dimensional  model.  The  results  of  these  tests  are 
presented  herein. 


PART  II:  TSUNAMI  CHARACTERISTICS  AT  CRESCENT  CITY,  CALIF. 

Maximum  Observed  Height  of  the  1964  Tsunami 

7.  Crescent  City  is  located  at  a  point  on  the  California  coast 
that  is  particularly  susceptible  to  the  damaging  effects  of  tsunamis.  The 
largest  tsunami  recorded  at  Crescent  City  was  the  one  generated  by  the 
1964  Alaskan  earthquake.  Fig.  1  is  a  reproduction  of  the  marigram  of  the 
1964  tsunami  taken  from  the  wave  gage  at  Citizen's  Dock  in  Crescent  City 
Harbor . 


Fig.  1.  Qagi5  record  of  1964  tsunami  at  Crescent  City 


8.  The  first  problem  was  that  of  organizing  the  information  con¬ 
tained  in  the  various  reports  concerning  the  observed  heights  of  the  1964 
tsunami  at  Crescent  City.  After  reviewing  the  available  data,  it  is  be¬ 
lieved  that  the  tide  gage  at  Citizen's  Dock  was  rendered  inoperable  by  the 
wave  crest  preceding  the  largest  wave.  Therefore,  in  order  to  assimilate 
the  tsunami  observations,  the  sequence  of  wave  crests  has  been  designated 


k 


as  shown  in  fig.  2.*  Accordingly,  tie  wave  crest  that  destroyed  the  tide 
gage  was  the  fourth  crest,  while  the  largest  wave  was  the  fifth  crest.  The 
dashed  lines  in  fig.  2  indicate  Wilson's  extrapolations  of  the  wave  record. 

9«  Since  it  was  considered  undesirable  to  rely  on  an  extrapolation 
of  the  crest  height  for  the  fifth  crest  from  the  existing  record,  obser¬ 
vations  of  eyewitnesses  were  the  only  available  means  of  estimating  the 
rise  and  fall  characteristics  of  the  fifth  and  largest  crest.  The  ob¬ 
servations  from  four  sources  are  tabulated  in  fig.  3*  It  is  interesting 
to  note  that  there  is  a  large  degree  of  disagreement  as  to  the  magnitudes 
of  the  crest  heights,  even^for  crests  that  were  recorded  by  the  tide  gage. 
This  may  be  due  in  part  to  the  disagreement  as  to  what  constituted 
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Fig.  2.  Wave  configura-  Fig.  3-  The  selected  design  wave 

tion  of  1964  tsunami 
(from  Wilson) 


*  Fig.  2  is  from  reference  1;  crest  notations  have  been  added. 

5 


✓ 


distinct  troughs  and  crests.  The  fifth  crest  and  preceding  trough  heights 

O  ii 

were  reported  by  WiegelJ  and  Magoon  after  extensive  interviews  with  re¬ 
liable  observers  who  were  on  or  near  Citizen's  Dock  during  the  tsunami. 

As  a  result,  the  reported  heights  should  be  within  a  tolerable  range  of 
accuracy.  For  this  report,  and  for  any  subsequent  model  testing,  an 
assumed  maximum  wave  with  the  height  characteristics  shown  at  the  bottom 
of  fig.  3  was  selected. 

Period  of  the  1964  Tsunami 

10.  Tudor, ^  in  reconstructing  the  Crescent  City  marigram  for  the 
1964  Alaskan  tsunami,  gives  the  following  times  (PST)  for  the  occurrences 
of  the  crests:  firsij.  crest  2352,  27  March;  second  crest  0020  and  fourth 
crest  0115,  28  March.  What  Tudor  calls  the  fourth  crest  has  been  indicated 
in  this  report  as  the  fifth  crest.  Accordingly,  the  perio^L  of  the  first 
wave  is  28  min  and  the  mean  period  of  the  sequence  of  the  first  four  waves 
is  27.6  min.  This  value  can  be  compared  with  the  periods  of  the  corre¬ 
sponding  tsunami  observations  on  the  northwestern  coast  of  the  United 
States.  These  periods  are  approximately  as  follows: 


Location 

Period,  min 

Location 

Period,  min 

Neah  Bay,  Wash. 

31 

Rincon  Island,  Calif. 

30 

Friday  Harbor,  Wash. 

31 

Santa  Monica,  Calif. 

26 

Astoria,  Or eg. 

30 

Los  Angeles,  Calif. 

30 

San  Francisco,  Calif. 

30 

Newport  Bay,  Calif. 

31 

Alameda,  Calif. 

30 

La  Jolla,  Calif. 

31 

Avila  Beach,  Calif. 

30 

San  Diego,  Calif. 

31 

The  above  values,  which  were  obtained  by  scanning  the  tsunami  tide  gage 

records  shown  by  Spaeth  and  Berkman,^  are  close  enough  to  the  periods  of 

waves  assigned  to  Crescent  City  so  that  the  value  of  T  =28  min  can  be 

P 

selected  as  the  characteristic  period  of  1964  Alaskan  tsunami.  Tide  gage 
records  of  the  1964  tsunami  at  other  localities  show  the  presence  of  waves 
of  other  periods  that  are  missing  in  the  corresponding  record  for 
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Crescent  City  owing  to  the  early  destruction  of  the  gage.  These  may  he  in¬ 
ferred  from  the  very  complete  gage  records  taken  at  Crescent  City  during 
the  i960  tsunami  from  Chile. 

Spectral  Analysis  of  a  Portion  of  the  i960  Tsunami  Records 

11.  A  useful  method  for  describing  a  vibrating  environment,  -whether 
it  be  random,  periodic,  or  a  mixture  of  both,  is  to  determine  the  spectral 
composition  of  the  response  data  representing  the  environment.  Several 
different  spectral  density  functions  might  be  conceived  and  employed  to 
describe  the  frequency  composition  of  random  or  periodic  vibration  re¬ 
sponses  of  -which  two  are  amplitude  spectra  and  power  spectra.  Each  of  the 
periodic  spectra  consists  of  discrete  -components  at  specific  frequencies 
and  harmonics  of  these  frequencies;  for  random  vibration  data,  each  of  the 
spectra  is  continuous ^with  responses  possible  at  any  or  all  frequencies  in 
the  spectrum. 

12.  Because  of  the  relatively  small  number  of  visually  obvious  dis¬ 
crete  frequencies,  the  wave  records  of  the  May  i960  tsunami  at  Crescent 
City  were  treated  as  periodic  vibration  data.  Any  periodic  signal  y(t) 
can  be  expressed  by  a  Fourier  series  as  follows:^ 

y(t)  =  Y0  +  Yi  sin  (Srrfjt  +  ftj)  +  Y2  sin  (Im^t  + 


+  Y_  sin  (6nf  t  +  ^ )  +  ... 


GO 

=  Y0  +  ZYn  sin  +  V 


where  is  the  fundamental  frequency  in  cps,  Y^  is  the  amplitude  of  the 
nth  component,  and  ^  is  the  phase  angle  of  the  nth  component.  A  de¬ 
scription  in  terms  of  the  harmonic  amplitudes  and  frequencies  is  possible 
by  frequency  spectrum  analysis  and  is  presented  as  the  remainder  of  this 
section. 

13.  Mr.  J.  L.  Ferguson  of  the  WES  Electronic  Data  Reduction  Lab 
undertook  the  task  of  producing  the  amplitude  spectra  (frequency  versus 
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wave-height  plots)  for  &  portion  of  the  Crescent  City  tsunami  height  record 
of  May  i960.  The  record  submitted  for  analysis  was  an  average  of  the  read¬ 
ings  taken  at  Citizen's  Dock  and  Dutton's  Dock  in  Crescent  City  Harbor. 

Ihe  average  record  was  done  by  hand  by  personnel  of  the  Hilo  Model  Branch, 
U.  S.  Army  Engineer  District,  Honolulu,  in  1965*  Because  of  the  2-min 
phase  lag  between  the  response  at  the  two  locations,  each  of  the  records 
was  offset  1  min  before  the  average  heights  were  computed.  The  results  of 
this  effort  are  presented  as  figs.  5*8. 

14.  The  si2ps  taken  in  processing  the  pictorial  wave  record  prior  to 
reduction  by  the  analog  wave  analyzer  are  outlined  below: 

a.  The  pictorial  record  of  the  tsunami  was  transferred  to  a 
~  digital  magnetic  tape  (in  a  compatible  BCD  format)  by  a 

pencil-line  follower  (Edwin  Industries  model  PF1QA).  The 
follower  digitized  and  recorded  the  y-coordinates  (frou  a 
prescribed  base  line)  at  1-nm  horizontal  increments  along 
the  entire  18-5 -hr  pictorial  record. 

b ,  A  digital  to  analog  conversion  was  then  effected  by  a  small 
digital  computer  (Scientific  Controls  Corporation  irodel 
650-2),  and  the  resulting  analog  signal  was  processed  by  tl.e 
wave  analyzer  (Gulton  Ortholog  model  OR-WA/l)  diagrammed  in 
fig*  4. 

15.  Each  of  the  frequency  spectra  (figs.  5-8)  was  generated  by  a 
slow-moving  filter  or  "window"  of  the  indicated  bandwidth.  The  entire 


•  Mteraerih  fitters  ae  diseased  in  mst  caaotchensive  texts  deilin*  wift  tto  ttwxy  d  electrics!  filters. 

Jig.  4.  Analog  wave  analyzer 
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Fig.  5.  Frequency  spectra  for  entire  tsunami  record,  0530  to  2400  hr 
23  May  i960,  at  Crescent  City  with  0. 00194- ,  O.OO388-,  and  0.0097- 

cpm  filters 
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Fig.  6.  Frequency  spectra  for  the  period  0530  to  1300  hr,  23  May  i960, 
at  Crescent  City  with  0.0019*+-,  0,00388-,  and  0.0097-cpra  filters 
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Fig.  7*  Frequency  spectra  for  the  period  1300  to  1700  hr,  23  May  i960, 
at  Crescent  City  with  0.0019*+- j  0.00388-,  and  0.0097-cpm  filters 
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Fig.  8.  Frequency  spectra  for  the  period  1700  to  2k00  hr,  23  May  i960, 
at  Crescent  City  with  0.0019^-,  O.OO388-,  and  0.0097-cpm  filters 
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analog  wave  record  signal  was  continuously  cycled  through  the  analyzer  and 
the  window  position  moved  so  slowly  across  the  spectrum  that  any  wave 
height  may  be  considered  as  being  associated  with  its  discrete  frequency. 
All  signal  amplitudes  of  a  particular  frequency  were  averaged  (by  means  of 
an  RC  averaging  network  within  the  wave  analyzer)  over  the  entire  time 
period  to  obtain  the  average  wave  height.  This  height,  as  shown  in  the 
plots,  would  not  be  comparable  to  a  maximum  height  as  measured  from  the 
wave  record.  Since  none  of  the  frequencies  appeared  at  all  times  on  the 
wave  record,  the  heights  shown  on  the  spectrum  plots  are  lower  than  the 
maximum  heights  measured  on  the  pictorial  record  with  the  heights  for  the 
least  occurring  frequencies  being  the  most  out  of  proportion. 

1 6.  Each  of  the  four  separate  record  lengths  of  the  i960  tsunami  at 
Crescent  City  was  subjected  to  analysis  using  three  separate  window  widths. 
Hie  entire  record  (from  0530  to  2k00  hr,  23  May  i960)  was  analyzed  first. 
The  record  was  then  divided  into  three  parts,  0530  to  1300  hr,  1300  to 
1700  hr,  and  1700  to  2k00  hr;  and  separate  analyses  were  performed  on  each 
of  these  parts.  Because  of  the  amplitude  averaging  technique  and  the 
distribution  of  the  wave  record  spectral  frequencies,  the  plotted  height 
increased  as  the  window  width  increased  for  the  three  separate  analyses  of 
the  four  series  of  runs. 

17.  The  statistical  methods  commonly  used  to  establish  confidence 
limits  with  respect  to  the  accuracy  of  -fee  wave-height  measurements,  shown 
in  figs.  5-8,  are  the  statistical  coefficient  of  variation  and  the  chi- 

Q 

square  distribution. 

18.  The  statistical  coefficient  of  variation  expresses  the  probable 

error,  e„  ,  as 
s 

e.  -  '  7==  (2) 

VBW  X  SL 

where 

eg  =  probable  error  inherent  in  the  estimate 
BW  =  effective  bandwidth  of  band-pass  filter 
SL  *  sample  length  (in  units  corresponding  to  the  bandwidth) 

Moody  states  that  this  method  should  not  he  used  when  the  numerical  value 
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20.  Confidence  levels  for  each  of  the  time  intervals  studied,  with 
their  corresponding  bandwidths,  are  shown  in  fig.  9*  A  large  number  of 
degrees  of  freedom  are  required  to  obtain  a  high  confidence  in  an  estimate. 
Equation  3  shows  that  this  can  be  accomplished  by  using  either  wider  filter 
bandwidths  or  longer  sample  lengths.  For  reasonable  sample  lengths  the 
problem  is  most  acute  at  the  low-frequency  end  of  the  spectrum  because  of 
the  narrowness  of  the  filter  required  to  resolve  the  resonances. 

21.  It  should  be  noted  from  fig.  9  that,  while  degrees  of  freedom  in 
excess  of  100  do  not  produce  proportionately  more  confidence,  those  obser¬ 
vations  with  less  than  20  degrees  of  freedom  produce  such  wide  confidence 
intervals  that  they  are  of  questionable  value. 

22.  While  one  cannot  place  any  confidence  in  the  fact  that  the 
plotted  heights  are  peak  values,  the  predominant  periods  of  80.0,  33*0, 
22.0,  and  18.0  min  are  clearly  defined  in  all  cases  and  can  be  accepted 
with  confidence.  !Ihese  predominant  periods  are  in  good  agreement  with 
those  mentioned  in  the  following  discussion  of  resonance,  which  were  com¬ 
puted  by  considering  the  oscillations  of  the  Crescent  City  area  from  the 
offshore  shelf  to  the  beach. 


Resonance 


23.  The  Omori  principle  that  the  tsunami  periods  observed  in  a  bay 
are  nearly  the  same  as  those  of  the  usual  secondary  oscillations  of  that 
bay  has  been  confirmed  by  Honda,  Terada,  and  Isitani.10  The  principle  is 
shown  to  be  valid  in  the  undulations  of  numerous  bays  along  the  coasts  of 
the  Pacific  Ocean  and  the  Sea  of  Japan.  Suda  and  Seki  in  their  study  of 
the  Atacam  tsunami  of  1922  and  of  the  Aleutian  tsunami  of  1923  found  that 
the  second  and  third  harmonics  of  the  fundamental  periods  are  highly  pre¬ 
dominant  in  bay  areas.11  Ihe  latter  finding  appears  to  apply  also  to  the 
undulations  of  the  i960  tsunami  at  Crescent  City. 

24.  The  Omori  principle  suggests  that  the  greater  tsunami  heights 
observed  in  some  bays  should  also  be  examined  on  the  basis  of  resonance. 
Wilson1  observes  that  the  excessive  heights  noted  at  Port  Alberni,  Canada, 
during  the  1964  Alaskan  tsunami  were  brought  about  because  of  the  parabolic 


variation  of  channel  depths ,  with  the  width  being  nearly  constant .  Like¬ 
wise,  the  greater  heights  observed  at  Port  Littleton,  Key  Zealand,  were  due 
to  the  uniformly  sloping  floor  of  the  channel.  For  the  analyses  of  these 
two  cases,  Wilson  makes  a  .judicial  application  of  Lash's  treatment  of  the 
oscillation  of  hays  of  such  configurations,  imposing  the  condition  that 
the  mouths  of  the  bays  are  nodal  lines. ^ 

25.  Similarly,  the  periods  of  oscillations  of  the  offshore  shelf  in 
the  proximity  of  Crescent  City,  as  shown  in  fig.  10,  have  been  considered. 
Hie  shelf  is  regarded  as  sloping  uniformly  with  the  depth  of  the  shelf 

Hc  *  350  ft,  and  the  length  of  the  shelf  L  *  1 6  nautical  riiles .  Assum¬ 
ing  that  the  depth  at  the  coast  is  nil,  computations  for  the  periods  of 
the  various  harmonics  give  values  of  89,  35,  18,  15,  and  13  min,  which 
curiously  enough  are  in  approximate  agreement  with  the  periods  found  in 
the  spectral  analysis  of  the  i960  tsunami.  Other  computations  were  also 
made,  assuming  that  the  shelf  at  the  coast  has  a  finite  depth.  However, 
the  effect  of  this  assumption  is  to  increase  the  value  of  the  fundamental 
periods.  A  more  thorough  treatment  of  shelf  oscillations  is  presented  in 
Appendix  D.  Noting  that  the  Shelf  environment  enclosing  Crescent  City  may 

he  regarded  as  approximating  a  sendelliptic  basin,  Wilson  makes  use  of  an 

12 

analysis  by  Goldsborough  and  finds  a  satisfactory  agreement  between  the 
observed  and  computed  periods  of  the  i960  tsunami. 

26.  Hilo  Bay  and  the  sea  vicinity  of  Crescent  City  are  both  suscep¬ 
tible  to  resonance.  Since  the  periods  of  the  fundamental  oscillations  of 
these  two  localities  are  different,  the  magnification  of  heights  in  each 
will  be  different  for  a  given  tsunami  of  a  definite  period.  Another  con¬ 
sideration  which  may  be  significant  is  the  following.  The  phenomenon  of 
resonance  is  most  conveniently  observed  with  standing  waves  in  rectangular 
basins.  Experience  shows  that  as  the  periods  of  the  forcing  oscillations 
are  decreased  from  an  initially  large  value  and  are  made  to  approach  the 
period  of  the  greatest  fundamental  oscillation  of  the  basin,  the  wave  mag¬ 
nification  in  the  basin  increases  steadily  to  a  limiting  value.  Then, 
when  the  period  of  the  forcing  oscillation  is  further  decreased  even  by  a 
small  value,  the  wave  magnification  decreases  markedly.  There  was  an 
evidence  of  this  in  the  Hilo  Bay  model  study.  As  the  main  object  of  that 


WATER  DEPTH  -  H,  FT 


(NAUTICAL.  MILES) 


DISTANCE  ALONG  MINOR  AXIS  -X  (NAUTICAL  MILES) 

Fig.  10.  Crescent  City  shelf  location  and  dimensions 
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14 


Fig.  11.  Wave  amplification  in  Hilo  Bay;  pilot  model  study 

study  was  to  examine  the  effect  of  model  distortion,  for  the  simplicity  and 
the  purity  of  observed  results,  the  bay  barrier  was  completed  forming  a 

\ 

truncated  bay.  The  average  results  of  the  magnification  ratio  —  , 

o 

h^  being  the  wave  height  at  the  barrier  and  ho  the  wave  height  at  the 

T 

bay  mouth,  from  all  the  tests  are  shown  in  fig.  11  as  a  function  of  ~- 

c 

(T  being  the  period  of  the  forcing  waves  and  the.  fundamental  period 

TO 

of  oscillation).  There  is  a  marked  peak  in  the  magnification  values  at 
T  .  Ar” 

—  =  0*99  •  Accordingly,  the  interval  ~  for  -which  the  waves  are  appre- 

C  AT  c 

ciably  magnified  is  about  =  0.075  •  In  the  same  figure  the  theoretical 

c  1 

magnification  curve  is  also  given. 

27 .  Eventually  the  question  will  be  asked  whether  there  is  a  dis¬ 
cernible  and  consistent  correlation  between  the  tsunami  heights  observed 
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at  Hilo  Bay  and  those  in  the  sea  approaches  of  Crescent  City.  The  above 
discussions  suggest  that  correlation  may  be  noted  for  the  situations  where 
the  tsunami  waves  have  periods  which  are  removed  from  the  periods  of  the 

fundamentals  of  the  localities  by  an  amount  of  about  AT  =  0.075T  or 
greater.  c 
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PART  III:  PREPARATION  OF  TSUNAMI  REFRACTION  DIAGRAMS 


Theoretical  Development 


28.  After  a  considerable  literature  search  it  was  deemed  desirable 
to  use  the  digital  computer  facilities  at  the  U.  S.  Any  Engineer  Waterways 
Experiment  Station  (WES)  to  prepare  the  refraction  diagrams.  In  January- 
1968,  the  main  computer  center  utilized  an  8000- space  main  memory  GE-225 
computer  and  an  off-line  Calcomp  750  incremental  plotter.  Attempts  were 
first  made  to  modify  existing  computer  programs  for  use  on  the  GE-225; 
however,  the  existing  programs  not  only  proved  difficult  to  modify  but  did 
not  incorporate  all  tie  features  considered  necessary  for  this  particular 
problem.  Therefore,  a  computer  program  was  written  specifically  for 
tsunami  refraction. 

29-  Ippen^  presents  the  basic  equation  for  wave  ray  refraction  as: 


dfl 

ds 


i  (sin  8  g-  cos  9g) 


(4) 


where,  considering  an  x  ,  y  rectangular  coordinate  system,  ^  represents 
the  change  in  the  angle  0  with  respect  to  the  distance  moved  along  the 
ray  (ds)  and  and  ^  represent  the  partial  derivatives  of  the  veloc¬ 
ity  in  the  x-  and  y-directions ,  respectively.  C  is  the  wave  celerity  and 
8  the  instantaneous  angular  orientation  of  the  ray. 

30.  Consider  two  positions,  Pq  and  P^  ,  on  a  ray  moving  in  the 
x-,  y-plan-a  as  shown  in  fig.  12;  PQ  and  P^  represent  two  successive 
points  on  the  ray  with  a  velocity  of  travel  at  each  point  of  C  =  i/gH  . 
Considering  H  as  the  depth  of  water  through  which  this  ray  travels, 


£  ,  I  4/2  H-l/?  Si 

&x  2  6  ax 

and 


i  M 

C  dx  2H  dx 


(5) 
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Hg.  12.  Wave  ray  notation 


In  a  like  manner, 


1  &  1_  Mi 
C  3y  2H  ay 

Substitution  of  equations  5  and  6  into  equation  h  yields 

sin  0I  -  cos  eS) 
ds  "  2H  \sln  9  3x  3y  / 

Noting  that  sin  0  =  ^  and  cos  0  =  ~  ,  equation  7  becomes 


d0  _  1/ 
ds  "  2H  l 


ds  ax  ds  by) 


Using  a  small  increment  for  ds  allows 


ds 


=  sj dx2  + 


dy 


(6) 


(7) 


(8) 


(9) 
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Substitution  of  equation  9  into  equation  8  yields 


d9  _  _  3H\  .  f _ 1 

ds  2H  vdx  dx  By/  / - 


LM«rj 

d0 

Small  shows  that  the  curvature  —  of  a  curve  y  =  f(x)  at  a  point 
P(x,y)  in  rectangular  coordinates  is  given  by: 


ds  =  r  ,13/2 

' 1  +  (£)J 

where  and  are  to  be  evaluated  at  P(x,y)  .  therefore, 

cbc 

equating  equation  10  to  equation  11  yields 


UA 

j- + (£f 


_ =  2H  _ L 

3/2  2H\dx3x  By/  , - 


■  (8  fi  -  f )  t  *  (»*]  ir 


A  Taylor's  expansion  of  the  function  of  (y  +  Ay)  yields 


dy  1  dy  2 

y  +  Ay  =  y  +TtAx  +  r  — £  ax  +  . . . 

o  o  dx  2,2 

dx 


or,  considering  only  the  first  three  terms, 


dx 

Substitution  of  equation  12  into  equation  13  yields 
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cjy  5h  3H 

Equation  l4  shows  that  by  starting  with  known  values  of  H  >  5  Bx  5  ciy 

at  Pq  ,  and  for  an  assigned  value  of  Ax  ,  one  can  compute  the  Ay  of 
point  P,  .  The  location  of  P1  is  then  simply  x  +  Ax  ,  y  +  Ay  . 

1  °  u  dy 

31.  Tc  move  from  ?1  to  point  V  ,  it  is  necessary  to  know  at 

P  .  go  SC  at  P,  is  tan  (0  +  A0)  and  it  can  be  shown  that 
1  9  dx  1  o 

tan  (0  +  A©)  =  tan  0Q  +  (1  +  tan2  0q)  a9  ^ ) 


Equation  10  can  be  written  as 


2H  V  dx  3x 


f * (if 


1  /d£  m  2H\  Av 
A0  =  2H  Vdx  Sx  9yj  ** 


Therefore 


or  from  equation  12 


(5x\  (5sc\  +  (lX) 

U/Pi  Wpo  w/p 


32.  The  time  of  travel  along  an  increment  As  is  given  by 


At  = 


(Ax2  +  AV2) 


Thus,  for  a  given  starting  position  (x,y),  initial  angle  ,  and  com¬ 
putational  increment  Ax  ,  the  corresponding  Ay  increment  can  be  calcula¬ 
ted.  This  next  point's  (PY  in  fig.  12)  coordinates  (x  +  Ax  ,  y  +  Ay) 
and  the  computed  then  become  the  basis  for  computation  to  P^  • 


This  process  is  continued  throughout  the  length  of  the  ray. 

33*  It  is  apparent  from  equation  14  that  it  is  necessary  to  deter¬ 
mine  the  depth  at  a  given  point  (such  as  P  )  in  order  to  progress  to  the 
next  point  (P-^) .  Equation  l4  shows  that  the  incremental  computation  of  the 
orthogonal  path  is  expressed  as  a  function  of  H(x,y)  ,  so  the  path  of  the 
ray  can  be  constructed  across  a  grid  of  water  depths,  each  having  an  x-  and 
y-dimension  in  the  grid  area. 

34.  A  regular,  square  grid  of  finite  boundaries  was  overlaid  on  a 
map  showing  the  area  through  which  the  waves  were  to  travel,  and  depth  val¬ 
ues  were  assigned  to  each  of  the  grid  intersections.  In  general,  however, 
the  points  for  which  a  depth  value  is  required  (beginning  points  of  each 
incremental  segment)  will  not  fall  on  the  regular  grid  points,  so  it  is 

necessary  to  interpolate  for  the  required  values  of  the  depth. 
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35  *  Dobson  utilized  a  least  squares  method  to  fit  a  surface  of 

12  grid  points  in  the  immediate  vicinity  of  the  area  of  interest  and  Wil- 
l8 

son  used  a  least  squares  fit  to  the  four  closest  depth  values  to  calcu¬ 
late  the  depth  at  an  arbitrary  point  within  a  grid  square.  It  is  apparent 
that  very  few  surfaces  can  be  accurately  described  by  a  second  degree  equa¬ 
tion  when  applied  to  the  entire  grid;  but  when  it  is  fitted  locally,  in  a 
mosaic  manner,  it  is  a  relatively  accurate  method  for  interpolating  the 
depth  at  intermediate  points.  Of  course,  the  method  would  yield  best  re¬ 
sults  if  the  ocean  bottom  was  a  second  order  surface,  but  should  yield 
satisfactory  results  if  a  small  enough  grid  is  used.  The  method  used  to 

^TT  air 

find  the  interpolated  depth  as  well  as  ~  and  is  outlined  in  the 

ox  ay 

following  paragra^Jis. 

36.  Fig.  13  shows  the  grid  points  surrounding  the  area  of  interest 
in  which  it  is  desired  to  find  the  depth  at  a  randomly  placed  point  P^  . 
This  can  be  done  if  the  depths  at  each  grid  intersection  are  known.  It  was 
decided  that  best  results  could  be  obtained  if  two  surfaces  were  fit 
through  the  area  of  interest  and  if  the  depth  were  computed  for  the  random  ■ 
point  P1  which  is  the  average  of  these  two  surface  fits;  likewise,  it 
was  considered  best  to  utilize  the  average  values  of  and  ~  computed 
from  each  of  these  two  surfaces. 


37.  One  may  solve  for  the  depth  at  point  P. 


i  ’  \ 


,  using 
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equation  19  where  the  coefficients  A  ,  B  ,  C  ,  D  ,  and  E  are  computed  by 
taking  the  average  values  of  two  sets  of  coefficients  obtained  from  two 
different  surface  fits  through  the  same  area. 

Hp  =  +  A  6x  +  B  6y  +  C  fix  6y  +  D  6x2  +  E  6y2  (19) 

38.  If  one  considers  the  depth  values  ^  ,  H^+2  ^  , 

Vl,3+D  ••  Vl  J-D  ’  “d  H(i,i.l)  as  6hCTm’ln  "«■  l4>  the  set’of 

simultaneous  equations  20  may  be  used  to  solve  for  the  first  set  of 
coefficients . 

39-  Another  surface  fit  through  the  area  of  interest  considering  the 
depths  shown  in  fig.  15  will  yield  another  set  of  coefficients  describing  a 
second  order  surface  passing  through  the  same  area  of  interest.  The  set  of 
simultaneous  equations  used  to  solve  for  a  second  set  of  coefficients  is 
Shown  as  equation  21. 


\  i«) 


H(»+i,  j+a) 
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Fig.  14.  Arrangement  of  grid  points  for  fitting  the  first 
surface  through  the  area  of  interest 


H(i+2,j)  =  H(i,j)  +2Alx  +  4V2 
H(i,j+1)  =  H(i,j)  +  +  Exy2 

H(i,j-1)  =  H(i,j)  -  \y  +  \y2 

H(i+i,j+D  =  H(i,j) +  v  +  V  +  cixy  +  v2  +  v2 

H(i-5-i,j-D  =  H(i}j) +  v  "  Biy  ■  cixy  +  Dix2  +  Eiy2 


£2t-A*EA  0F '"TEREST 


Fig.  15.  Arrangement  of  grid  points  for  fitting  the  second 
surface  through  the  area  of  interest 


>i,j)  =  H(i,j)  +  V +  V 


(i.j+D  =  Ha,j)  *V'V 
(i,J*2)  =  H(1,J)  *a&*  ‘‘V 


(i+2, j+1)  *  Hu,j)  +  M2x  +  V  +  *2*  *  4V2  *  VS 


(i+i, j+2)  -  H(i,j) ♦  V +  +  2v2  *  v2 
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40.  The  average  value  for  each  of  these  two  sets  of  coefficients  Is 
given  by 


A  = 


A1  +  AS 


B  = 


C  = 


D  = 


E  = 


31  +  B2 


C1  +  C2 


Dl  +  D2 


E,  +  E„ 

±  C. 


(22) 


J 


Substitution  of  these  coefficient  values  and  the  correct  values  for  the  lo¬ 
cations  (6x,6y)  of  the  randomly  spaced  point  in  the  grid  square  into  equa¬ 
tion  19  will  yield  the  depth  (Hp^ )  at  the  random  point.  Partial,  differen¬ 
tiation  of  equation  19  with  respect  to  x  end  y  ,  respectively,  will 
yield 


H  =  A  +  Cty 


?D6x 


and  (23) 

=  B  +  C6x  +  2E6y 

Computer  Program 

4l.  A  listing  of  the  program  with  an  explanation  of  the  program  var¬ 
iables  and  an  exauqple  of  input  data  are  presented  in  Appendix  A.  The  depth 
values  used  as  input  data  for  the  program  were  a  weighted  average  of  three 
independently  obtained  depths.  One  of  the  sources  used  was  the  ocean 

depths  averaged  over  areas  of  one-degree  squares  of  latitude  and  longi- 
19 

tude  and  compiled  ns  computer  printout  available  from  the  National  Ocean¬ 
ographic  Data  Center,  The  second  and  third  sources  of  depth  data  (both 
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from  the  U.  S.  Naval  Oceanographic  Office)  were  an  unpublished  collection 
of  ocean  sounding  sheets  at  an  approximate  scale  of  1:800,000,  and  the 
H.O.  series  of  charts  at  an  approximate  scale  of  1:5,990,000.  The  plots 
that  were  obtained  from  the  computer  were  referenced  to  the  1:5, 990*00° 
series  charts. 

42.  The  presentation  of  a  geoidal  surface  (the  earth)  on  a  map  re¬ 
quires  that  the  representation  be  distorted  in  some  way.  Mosi  Mercator 
projections  (the 
most  common  chart¬ 
ing  method  for  maps 
used  for  navigation  „  so 

and  ocean  studies)  « 

'  O 

IU 

are  projected  in  40 

true  scale  at  0°  ° 

t 

latitude  (the  equa-  £  so 

J 

tor)  with  the  dis-  * 

K 

tort ion  becoming  °  20 

greater  the  farther 

10 

one  moves  from  the 
equator  as  shown  in  ° 

fig.  16.  At  a 

10 

given  point  on  a 

Mercator  map  the  3  20 

111 

distortion  is  the  o 

bJ 

same  in  all  direc-  30 

O 

tions  and  is  equal  ? 

to  the  secant  of  <  40 

the  latitude  at  jj 

.0 

which  one  is  lo-  '  “  » 

cated  (i.e.,  at  60° 

north  or  south 

•  60 
latitude  a  Mercator  140  150  160  170  100  170  160  150  uo  130 

.  . .  ,  .  .  .  EAST  LONGITUDE,  DEGREES  WEST  LONGITUDE,  DEGREES 

map  is  distorted  by 

a  factor  of  2  0)  Fig*  F<iual  distances  represented  at  different 


IU  20 
in 


140  150  160  170  160 

EAST  LONGITUDE,  DEGREES 


.70  160  150  140  «30 

WEST  LONGITUDE,  DEGREES 


Fig.  l6.  Equal  distances  represented  at  different 
latitudes  on  a  Mercator  projection 
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To  get  a  true  picture  of  wave  travel  over  large  distances  and  particularly 
through  many  degrees  of  latitude,  one  must  compensate  for  this  distortion. 

43.  Two  methods  of  compensation  were  considered  for  this  study. 

Dr.  Ledolph  Baer  of  the  Lockheed-Califoraia  Company  has  devised  an  essen¬ 
tially  undistorted  flat  representation  of  the  earth's  surface  by  drawing 
a  map  on  an  icosahedral-gnomonic  projection.  The  earth's  surface  was 
divided  into  20  sections  and  each  of  these  sections  was  individually  pro¬ 
jected  (using  a  gnomonic  projection)  onto  a  flat  surface  with  the  center 
point  of  each  section  being  the  tangent  point  between  the  earth  and  the 
map.  Each  of  the  sections  is  therefore  undistorted  at  its  center,  but  be¬ 
comes  increasingly  distorted  toward  the  edges  of  that  individual  section. 

44.  The  second  method  studied,  and  the  one  adopted  for  application 
to  this  problem,'  was  to  distort  each  computed  or  assigned  A  distance  by 
the  secant  of  the  latitude  of  its  location.  Careful  examination  of  the 
theoretical  treatment  of  the  problem  revealed  that  this  could  be  easily 
accomplished  by  multiplying  each  of  the  depth  values  which  were  to  be  input 
to  the  program  by  the  square  of  the  secant  of  the  latitude  at  which  these 
depth  values  were  located.  A  rigorous  treatment  of  this  is  presented  in 
Appendix  B.  This  multiplication  was  performed  before  inputting  the  depth 
data;  therefore,  all  the  depth  values  of  the  input  data  grid  shown  in  Ap¬ 
pendix  A  have  already  been  multiplied  by  the  secant  squared  of  the  latitude 
of  their  respective  location.  This  has  the  effect  of  showing  the  wave  path 
on  the  map  as  if  it  had  taken  place  on  the  earth's  surface  and  had  then 
been  transferred  point  by  point  to  the  map. 

45.  To  test  the  foregoing  numerical  techniques,  it  was  deemed  neces¬ 
sary  to  compare  the  computer -dr awn  refraction  diagrams  with  those  obtained 
by  a  closed  solution  of  wave  refraction  around  a  given  island  shape.  Such 
an  example  is  described  in  Appendix  C. 

Refraction  Diagrams 

46.  The  ref /action  diagrams  computed  and  drawn  by  the  previously  de¬ 
scribed  technique  are  shown  in  figs.  17-25.  The  figures  show  three  recent 
tsunami -front  patterns  (Alaska,  1964;  Chile,  i960;  and  Kamchatka,  1952) 
from  each  of  the  given  sources  to  Crescent  City  and  Hilo,  Hawaii.  The 
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40' 


I25°00 


Fig.  21.  Refraction  diagram  in  Crescent  City  area  from  Chilean  tsunami  of  May  22,  i960 


Fig.  22.  Refraction  diagram  in  Hilo  area 
from  Chilean  tsunami  of  May  22,  i960 
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Fig.  23.  Refraction  diagram  for  Kamchatka  tsunami  of  November  4,  1952 


125°00 


Fig.  24.  Refraction  diagram  in  Crescent  City  area  from 
Kamchatka  tsunami  of  November  4,  1952 


waves  were  drawn  threap  the  open  ocean  areas  cm  1:6,000,000- scale  naps  and 
upon  arrival  in  the  nearshore  areas  of  Ctes cent  City  and  Hilo,  the  tsunami 
fronts  \«sre  transferred  to  a  larger  scale  map  (l:196,S48  in  the  case  of 
Crescent  City  and  1:1,650.000  in  the  case  of  Hilo)  and  continued  to  very 
near  the  two  respective  points  of  interest. 

4?  *fhe  tsunami  fronts  from  each  of  the  three  locations  are  clearly 
defined  in  the  Crescent  City  area,  thus  describing  definite  wave  fronts 
for  model  wave  generator  plactfnent,  Hie  tsunami -front  patterns  approach 
co-'nnidence  from  the  three  sources  as  they  approach  Crescent  City. 

48.  Refraction  theory  dictates  that  the  wave  height  h^  (shown  in 
fig.  26)  will  be 


fiY 

\\/ 


where 

hQ  =  known  wave  height  at  any  distance  x  from  the  epicenter 

bc  =  physical  distance  between  rays  at  some  known  wave  height  (hQ) 

b,  =  physical  distance  between  rays  at  a  point  at  which  one  is 
seeking  a  wave  height 

Assuming  that  hQ  is  a  constant  known  quantity  for  each  of  the  beginning 
tsunamis  at  a  given  location,  the  quantity  (Wt^)1/2  should  give  the  rel¬ 
ative  wave  heights  at  any  other  desired  location  for  the  same  tsunami. 

The  distances  t)  and  b,  were  scaled  from  the  small-scale  sheets,  and 

0  ,i/2  b 

the  quantity  (b^/h^)  '  for  the  nearshore  areas  at  Hilo  and  Crescent  City 
is  shown  in  the  following  tabulation: 


Tsunami 

Alaska 

1964 


Hilo,  Hawaii 
/b  Y1/2 
(\)  =  °'5 


Crescent  City,  Calif. 


/b  X1/2 

(b2)  =  U* 

W 


Chile 

0.2 

0.09 

o 

& 

r-i 

Kamchatka 

0.6 

0.08 

1952 

It  must  be  noted  that  one  could  obtain  different  values  for  the  ratio 
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Fig.  2 6.  Refraction  notation 

(bo/bb)1/2  by  choosing  another  point  such  as  halfway  between  the  tsunami 

source  and  the  point  of  interest  for  the  measurement  of  h  and  b  . 

o  o 

Since  the  tsunami  height  hQ  probably  does  not  vary  much  in  the  open 
ocean,  the  selection  of  any  other  point  will  undoubtedly  change  the  final 
nearshore  wave  height  h^  by  changing  the  (^/b^)1/2  ratio. 

49.  No  attempt  was  made  to  present  a  design  wave  height  .  The 
four  lowest  values  of  (b^j/b^)1^2  given  in  the  above  tabulation  dictate 
that  the  wave  was  exceedingly  small  (considering  its  small  height  in  the 
open  ocean),  a  condition  that  was  not  present  at  either  of  the  two  given 
locations  for  any  of  these  tsunamis.  The  magnitude  of  these  terms  does 
infer,  however,  that  no  wave  magnification  takes  place  in  the  open  ocean 


due  to  gross  wave  refraction.  Local  severe  discontinuities  will  undoubtedly 
cause  local  wave  amplitude  increases  even  in  the  open  ocean;  however,  an  un¬ 
usually  large  amplitude  associated  with  any  tsunami  is  probably  due  to 
severe  shore  refraction,  shoaling  conditions,  and  the  resonance  of  the  par¬ 
ticular  nearshore  area.  Comparison  of  tsunami  heights  at  various  loca¬ 
tions  mentioned  in  Part  IV  shows  wave  heights  for  which  all  wave-height  mag¬ 
nification  factors  (refraction,  shoaling,  and  friction)  have  already  been 
considered.  It  is  seen  that  tsunamis  of  equal  periods  passing  through  equal 
water  depths  will  exhibit  equal  wave- height  increases  due  to  shoaling.  The 

time  of  travel  of  the  tsunamis  shown  in  figs.  17" 25  is  very  close  to  that 

6  ?l  PP 

presented  in  the  literature.  5  * 


PART  IV:  PROBABILITY  OP  OCCURRENCE  OF  TSUNAMIS  OF 
VARIOUS  HEIGHTS  AT  CRESCENT  CITY 


■Distribution  of  Tsunamis  of  Various 
Heights  at  Crescent  City 


50.  The  relation  of  the  tsunami  height  versus  frequency  of  occur¬ 
rence  for  Hilo  has  been  determined  by  Professor  Cox,  both  on  the  basis  of 
visual  observations  and  marigraphic  tide  recordings.  The  first  group 
covers  a  period  of  about  130  years  and  the  latter  a  shorter  period  of 
about  24  years.  The  heights  refer  to  runup  values,  which  are  derived  from 
the  given  bay  heights  by  multiplying  the  latter  values  by  1.5.  A  signifi¬ 
cant  fact  recognized  by  Cox  is  that  the  frequency  relation  needs  to  be  ex¬ 
pressed  differently  depending  on  the  magnitude,  of  the  heights.  For  the 
smaller  heights  a  power  law  is  indicated,  whereas  for  the  larger  heights 


the  logarithmic  form 
is  more  appropriate. 
Our  present  Interest 
is  in  the  latter. 
Keeping  in  mind  the 
frequency  relation 
for  bay  tsunami 
heights,  the  Cox  de¬ 
terminations  were  di¬ 
vided  by  1.5  and  are 
presented  in  fig.  27. 
Wiegel  gives  some 
similar  results  for 
Crescent  City  based 
on  16  observations-3 
of  which  those  relat¬ 
ing  to  larger  heights 
are  also  plotted  in 
fig.  27.  The  curve 
representing  average 
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Fig.  27.  Distribution  function 
for  maximum  wave  height 
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that  a  tsunami  will  strike  Crescent  City  in'  a  100-year  interval  with  a  max¬ 
imum  wave  height  of  at  least  23  ft;  furthermore,  there  is  a  20  percent 
chance  that  such  a  wave  will  occur  in  a  45-year  period,  and  a  95  percent 
chance  that  such  a  wave  will  strike  the  town  within  a  period  of  6 00  years. 
Agreeing  with  Wiegel,  it  will  be  remarked  that  an  extension  of  tsunami 
heights  to  values  much  above  those  that  have  been  observed  may  be 
erroneous . 


Coaparison  of  Tsunami  Heights  at  Hilo  and  Crescent  City 


52.  In  the  treatment  shown  above,  a  tacit  assumption  is  made  that  a 
correlation  exists  between  the  tsunami  heights  at  Hilo  and  Crescent  City 
for  a  given  instance  of  tsunami  occurrence  ami  that  the  correlation  is  one 
of  equality.  Table  2  shows  the  heights  observed  at  Hilo  and  the  heights 
observed  at  Crescent  City  or  values  reduced  from  the  heights  observed  at 
San  Francisco  and  the 
neighboring  locali¬ 
ties  of  Los  Angeles 
and  San  Diego.  On 
the  basis  of  data  ob¬ 
tained  since  1948,  it 
is  found  that,  on  the  t 

O 

average,  the  Crescent  > 

6 

City  tsunami  heights 
are  2.7  timeB  greater  5 

M 

than  the  values  ob¬ 
served  in  San  Fran¬ 
cisco,  Los  Angeles, 
and  San  Diego  (see 
fig.  29).  The  source 
of  the  data  in 
table  2  relating  to 
tsunamis  1  to  13,  in¬ 
clusive,  is  a  catalog 


Fig.  29,  Correlation  of  tsunami  heights  between 
San  Francisco  and  Crescent  City 
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prepared  for  the  International  Tsunami  Information  Center*  Additional 

sources  of  data  for  tsunamis  l4  to  18,  inclusive,  are  reports  of  the  U.  S. 

6  21  22 

Coast  and  Geodetic  Survey.  *  *  Hie  heights  refer  to  the  maximum  rise 

or  fall  noted  during  a  tsunami.  The  values  from  table  2  are  plotted  in 
fig*  30.  If  the  points  relating  to  tsunamis  2,  3,  4,  and  5  in  table  2  are 
ignored,  the  remaining  points  suggest  that,  on  the  average,  during  a  given 
tsunami  the  wave  heights  at  Hilo  and  at  Crescent  City  are  similar. 


Fig.  30.  Respective  tsunami  heights  at  Crescent  City  and 
Hilo  for  the  same  tsunami 
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Irregular  Tsunami  Heights  at  Hilo 


53.  It  would  not  be  amiss  to  reconsider  tsunami  cases  2,  3,  4,  and 
5,  which  show  excessively  large  values.  These  are  runup  heights  and  even 
if  an  allowance  is  made  for  tWs  fact,  the  reduced  bay  height  values  are 
still  far  higher  than  might  be  expected.  Prom  a  tabulation  given  by 
Spaeth  and  Berkman,^  it  is  inferred  that  the  ratios  of  the  tsunami  heights 
at  Hilo  to  those  at  Honolulu  for  the  incidences  of  1952,  1957>  i960,  and 
1964  are  1.79>  2.78,  1.75>  and  4.62,  respectively.  The  mean  ratio,  2-73> 
may  be  used  to  infer  the  probable  tsunami  heights  at  Hilo.  This  reduction 
is  done  in  table  3  for  tsunamis  2,  3,  4,  and  5  (as  listed  in  table  2),  and 
it  is  seen  that  the  heights  inferred  for  Hilo  are  of  the  same  order  of 
magnitude  as  the  heights  at  Crescent  City.  Occasionally  Hilo  shows  exces¬ 
sively  high  runup  values  for  reasons  difficult  to  assess. 

Tsunami  Height  Attenuation  with  Distance  Traveled 

54.  The  data  points  of  fig.  30  (a  comparison  of  tsunami  heights  at 
Hilo  and  Crescent  City)  show  scattering.  This  possibly  could  be  explained 
on  the  basis  of  the  distances  of  the  observation  localities  from  the  tsu¬ 
nami  generating  area.  Hilo  is  nearly  a  central  poi at  of  the  Pacific  perim¬ 
eter  of  earthquake  activity.  It  is  thus  proper  that  the  tsunami  heights  of 
any  locality  be  compared  with  those  at  Hilo.  For  other  points,  the  effect 
of  distances  of  tsunami  travel  needs  to  be  considered.  This  is  particu¬ 
larly  true  for  localities  that  exhibit  resonant  conditions.  Distance 
affects  the  periods  which,  in  turn,  modify  the  magnification  of  waves  in 
bays  subject  to  resonance.  Distance  also  modifies  the  tsunami  heights. 

This  matter  is  examined  further  in  the  next  section. 

55*  In  order  to  accurately  correlate  wave  heights  at  Hilo  and  Cres¬ 
cent  City,  it  is  necessary  to  determine  the  dependence  of  wave  height  on 
travel  distance.  Unfortunately,  this  functional  dependence  has  been  deter¬ 
mined  for  only  very  ideal  cases  and  even  there  a  certain  amount  of  con¬ 
fusion  exists. 

56.  The  attenuation  of  wave  height  follows  the  general  relation 
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where 


h^  *  wave  height  at  a  distance  x  from  the  source 

h  =  wave  height  at  the  source 
o 

The  relation  is  unmanageable  in  this  form  since  the  wave  does  not  travel  in 
a  straight  line  from  point  to  point.  The  above  equation  was  derived  from 
measurements  and  theory  for  constant  average  depth  conditions  and  the  ap¬ 
proximate  values  for  the  exponent  n  are  tabulated  below. ^ 


Case 

n 

Case 

n 

Three-dimensional : 

5/6 

Two-dimensional : 

1/3 

Leading  waves 

Leading  waves 

Wave  body 

1 

Wave  body 

1/2 

Since  the  tsunami -generating  earthquakes  occur  along  a  fault  line,  the  wave 
form  may  actually  be  somewhere  between  these  two  extremes. 

57-  Tsunami  wavelengths  are  of  such  magnitude  that  tsunamis  may  be 
considered  shallow-water  waves  and  the  amplitude- to-depth  ratio  is  small 
enough  to  be  negligible.  The  wave  speed  is  thus  very  nearly  equal  to  the 
limiting  value 

C  *  (28) 

where 

C  ■  wave  celerity 
g  *  the  gravitational  constant 
Hq  »  the  depth  of  the  water 

One  would  expect,  therefore,  to  obtain  a  degree  of  correlation  between  wave 
height  and  time  of  travel  shown  by  the  general  form 

\  -  hx  t~°  (29) 

where,  in  this  case,  h^  is  the  wave  height  after  one  hour  of  travel,  t  is 
the  travel  time  in  hours,  a  is  a  curve  fitting  coefficient,  and  h^  is 
the  wave  height  after  t  hours  of  travel.  Available  data  for  the  i960 


tsunami  which  originated  off  the  Chilean  coast  are  plotted  in  figs.  31- 34 . 
The  travel  time  on  all  the  figures  was  determined  by  the  time  of  arrival 
of  the  first  wave  as  detect  d  by  tidal  gaging  stations.  Attempts  were  made 
to  correlate  time  with  both  che  initial  wave  height  and  the  maximum  oscil¬ 
lation  of  the  water  surface.  Figs.  31  and  33  exhibit  a  very  high  degree  of 
scatter,  so  the  values  were  averaged  over  convenient  intervals  and  plotted 
on  log- log  paper.  A  straight  line  was  fitted  to  the  log-log  plot  by  the 
least  squares  method  in  the  general  form  of  equation  29  and  the  resulting 
plots  are  shown  as  figs.  32  and  34. 

58.  The  wave  height  is  seen  to  diminish  with  an  increase  in  time 
according  to  the  general  exponential  relation;  however,  the  data  are  too 
erratic  to  use  as  a  basis  for  predicting  wave-height  distribution.  It  is 
generally  believed  that  the  nearshore  conditions  have  a  much  greater  effect 
on  the  wave  height  than  travel  distance.  Wave  measurements  would  have  to 
be  made  in  the  open,  ocean  to  avoid  the  complicating  factors  due  to  near¬ 
shore  and  harbor  topography. 


59 •  In  fig.  30  Crescent  City  wave  heights  are  compared  with  those 
of  Hilo  Bay,  and  the  Crescent  City  height  from  the  1964  Alaskan  tsunami 
is  close  to  the  first  rise  of  the  water  surface.  Importantly,  the  largest 
wave  height  observed  was  that  of  the  fifth  wave;  it  rose  about  23  ft,  and 
indicates  a  very  large  disparity  from  the  curve  purportedly  showing  the 
trend  in  the  correlation  of  heights  between  Hilo  Bay  and  Crescent  City. 
Uhls  large  disparity  cannot  be  explained  on  the  basis  of  travel  time  from 
the  generating  area  to  these  localities  under  consideration.  If  one  con¬ 
siders  hQ  and  to  be  the  tsunami  heights  at  the  shelf  approaches  to 

Crescent  City  and  to  Hilo  Bay,  respectively,  and  T  and  T.  to  be  the 

c  n 

travel  times  of  the  tsunami  wave  from  the  generating  area  to  these  locali¬ 
ties,  then  according  to  the  relation  shown  in  the  previous  section, 
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Fig.  32.  Log- log  plot  of  information  in  fig.  31 


3^.  Log-log  plot  of  information  in  fig.  33 


Spaeth  and  Berkman^  give  =  325  min  and  Tc  *  243  min,  and  hence 
hc/hfa  =  1.27  •  Refraction  analysis  shears  that  the  energy  reductions  at 
these  two  localities  are  also  of  approximately  similar  values.  Ihus,  one 
needs  to  seek  other  factors  for  the  anomaly.  One  possibility >  of  course, 
could  be  the  presence  of  disparities  in  the  resonant  stages  of  the  two 
localities  during  this  particular  tsunami.  One  would  expect  a  higher  wave 
magnification  at  Crescent  City  than  that  at  Hilo  Bay,  if  at  Crescent  City 
T  -  T 

the  quantity  — ^ — —  was  almost  nil,  and  at  Hilo  Bay,  of  a  moderate 
cr 

amount.  As  precise  determinations  of  the  critical  periods  of  the  two 

localities  are  as  yet  not  forthcoming,  the  matter  cannot  now  be  resolved. 
Another  possibility  could  be  an  uneven  concentration  of  energy  along  the 
crest  of  the  wave  front  at  the  generating  area.  This  is  also  a  matter 
that  is  not  amenable  to  an  accurate  determination. 

60.  Wilson  has  summarized  much  useful  information  relating  the 
fault  length,  tsunami  source  dimension,  tsunami  height,  and  tsunami  period 
and  frequency  with  earthquake  intensity.1  Such  information  gives  insight 
into  the  mechanism  of  tsunamis,  but  this  is  hardly  adequate  to  predict  what 
the  tsunami  heights  at  Crescent  City  would  be  if  another  earthquake  oc¬ 
curred  in  the  Alaskan  area  having  an  intensity  of  8.5  on  the  Richter  scale. 
The  most  relevant,  information  desired  is  the  frequency-height  relation  for 
Crescent  City.  From  the  latter  point  of  view,  the  exaggerated  disturbances 
of  the  1964  Alaskan  tsunami  hardly  fit  the  course  of  expected  events  of 
Crescent  City,  and  one  is  forced  to  regard  the  manifestation  as  a  singular 
event  of  rare  probability. 
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PART  V:  FUJMfi  TESTS  AND  RESULTS 

6l.  Because  of  the  long  period  and  large  height  of  the  1964  Crescent 
City  tsunami,  it  was  deemed  necessary  to  make  flume  studj.es  before  final 
selection  of  the  Crescent  City  model  configuration. 

Wave  Generator  and  Flume 


62.  An  existing  5-ft-wide,  100-ft-long  concrete  flume  was  used  for 
preliminary  testing.  The  object  of  the  tests  was  to  determine  the  charac¬ 
teristics  of  a  long-period  wave  model. 

63.  A  sketch  of  the  pneumatic  wave  generator  used  in  the  majority 
of  flume  tests  is  shown  in  fig.  35-  It  consists  of  a  rectangular  steel 
vessel  with  a  nozzle  attached  to  the  mouth  of  the  tank.  The  top  is  pro¬ 


vided  with  a  valve  lead¬ 
ing  to  the  aspirator  and 
another  set  of  valves  lead¬ 
ing  to  the  open  air.  In 
a  simple  manner  of  operation, 
the  flume  water  is  raised 

from  an  initial  level  H  , 

o 

measured  from  the  flume 
bottom,  to  a  level  H  . 
During  this  process  the  air 
valves  are  closed.  Next, 
the  aspirator  valve  is 
closed  and  the  air  valves 
are  opened.  Assume  that 
the  initial  wave  of  the 
efflux  is  rectangular  in 
shape  and  is  of  constant 
height  h  .  In  the  tenta¬ 
tive  design  of  the  tank,  use 
was  made  of  the  relation2^ 


3"  NIPPLE  CONNECTED  TO  VACUUM 
2"  SOLENOID  VALVE  — x 
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Fig.  35.  Pneumatic  wave  generator  used 
in  the  majority  of  flume  tests 


which  contains  in  addition  to  Hq  ,  h  ,  and  H  ,  the  quantities  k^  and 
X  .  The  significance  of  k^  is  shown  in  the  following  statements.  Prior 
to  the  opening  of  the  air  valves  the  suction  pressure  in  the  tank  is  &pQ  , 
say 

Apo  =  pg  (H  -  Ho)  (32) 


Following  the  opening  of  the  air  valves,  the  suction  pressure  is  decreased 
to  ,  allowing  the  flow  of  water  from  the  generator  and 


(33) 


The  value  of  k.,  may  vary  between  0.06  and  0.10  by  a  proper  selection  of 

the  air  gap  in  the  tank  and  the  aperture  sizes  of  the  valves.  A  pood  value 

to  take  is  k,  =  0.085  •  The  quantity  X  is  related  to  the  energy  loss  in 

26 

the  area  of  the  nozzle  and  would  have  a  reasonable  value  of  0.4. 

64.  The  generator  was  designed  with  the  present  1:125  undistorted 
scale  model  of  the  Crescent  City  coast  and  the  sea  environment  in  mind. 

The  prototype  depth  of  water  in  the  approaches  to  the  harbor  is  approxi¬ 
mately  62  ft.  This  gives  a  model  depth  Hq  of  0. 5  ft.  A  prototype  design 
wave  of  30  ft  would  have  a  value  h  =  0.24  ft.  Entering  equation  31  with 
the  quantities  h  =  0.24  ft,  Kq  =  0.5  ft,  X  =  0.4,  and  k^  =  O.O85  yields 
H  =  5.0  ft.  Allowing  1  ft  for  the  air  space,  the  height  of  the  generator 
tank  would  be  6  ft.  To  determine  the  length  of  the  tank,  it  is  imagined 
that  after  the  efflux  from  the  tank  is  completed,  the  generated  wave  is 
triangular  in  shape  and  extends  from  the  generator  to  the  model  coast. 

This  distance  is  about  80  ft  in  the  model.  At  the  front,  the  height  of 
wave  is  0.5  ft.  Thus,  the  wave  area  is  20  sq  ft,  requiring  that  the  length 
of  the  tank  be  4  ft. 


Model  Scale  Relations 


65.  The  derivation  of  scale  relations  becomes  more  meaningful  if 
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r  * 


Thus,  if  the  model  he  constructed  affinely  similar  to  the  prototype,  that 
is  at  corresponding  points 


6=6 
m  p 


then,  at  corresponding  points 


u  =  u 
m  p 


and 


provided  that 


which  is  the  scale  for  the  periods.  Let  S  be  the  scale  for  depths  and 
r  the  model  distortion,  a  quantity  larger  than  unity.  Then,  with 


S  = 


1 


rS  =  =£ 
m 


one  has 


66.  Three  series  of  tests  were  made  for  this  part  of  the  study.  The 
first  series  was  made  using  a  horizontal  and  vertical  scale  of  1:125.  This 
is  the  scale  of  the  existing  Crescent  City  model  which  was  constructed  at 
the  WES  for  short-period  wave  studies.  The  horizontal  scale  was  then  en¬ 
larged  to  1:375  (the  vertical  scale  remaining  1:125)  and  the  contours  in 
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the  flume  floor  remolded  for  the  second  series  of  tests.  A  set  of  hardware- 
cloth  wave  screens  was  then  placed  in  front  of  the  wave  generator  and  the 
third  series  of  tests  was  run  using  the  same  model  scales  as  those  used  in 
the  second  series.  5fce  flume  bottom  configurations  for  the  different  series 
of  tests  are  shown  in  fig.  3 6. 


O  6,000  10.000 


PROTOTYPE  DISTANCE,  UN  DISTORT  ED  FT 


Fig.  36.  Flume  bottom  configurations  for  the  three  series  of  tests 


First  Series  of  Utests,  Uhdistorted  Model 

67 .  A.  typical  wave  record  for  the  case  of  an  undistorted  model  with 
a  scale  of  1:125  is  shown  in  fig.  37.  This  series  of  tests  was  performed 
to  investigate  the  possibility  of  utilizing  the  Crescent  City  model  pres¬ 
ently  situated  at  the  WES.  It  became  evident  frcm  this  series  of  tests 
that  it  was  not  possible  to  develop  a  sufficiently  long-period  model  wave 
due  to  the  Interference  of  the  resonant  waves  in  the  flume.  The  resonant 
characteristics  of  this  size  model  cause  undesirable  interaction  which 
makes  it  impossible  to  model  even  a  half-period  wave. 


Second  SerieB  of  Tests .  Distorted  Model  Without  Screens 


68.  A  typical  wave  record  for  this  case  is  shown  in  fig.  38.  There 
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Fig*  37*  Typical 
wave  record  in 
undistorted  flume 
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Fig.  38.  Typical 

wave  record  in 

distorted  flume; 

T  =  35  sec 
m  J 


WAVE  ROD  AT  APPROXIMATE  LOCATION  OP  CITIZEN'S  DOCK 


^  WAVE  ROO  LOCATEO  AT  GENERATOR  fcvOL'TM 


t«  SO  SO  40  |0  SO 

MODEL  TIME,  SEC 

s  *  •«  FT  »s.n 

*«0  L  *  SS>S  FT  r  .  ss  SEC 


appear  to  be  two  objectionable  characteristics  in  this  case. 

a.  Bie  wave  record  is  very  spiked  with  the  initial  stage  of  the 
wave  being  greatly  affected  by  undesirable  exit  conditions 
of  the  wave  generator. 

b.  The  resonant  characteristics  of  the  model  again  make  it  im¬ 
possible  to  input  a  wave  of  the  desired  period.  The  flume 
configuration  consisted  of  a  long  uniform  section  up  to  the 
contoured  area  at  the  end  of  the  flume.  A  "coupling  effect" 
was  evident  with  the  two  areas  resonating  independently  at 
their  own  characteristic  frequencies  rather  than  as  a  whole. 


Third  Series  of  Tests,  Distorted  Model  with  Screens 


69.  Hardware- cloth  wave  screens  were  installed  in  the  distorted 
model  in  an  attempt  to  eliminate  the  objectionable  peaks  mentioned  in  the 
second  series  of  tests.  A  typical  wave  record  presented  as  fig.  39  shows 
that  this  objective  was  partially  accomplished;  however,  the  resonant  in¬ 
terference  the  barrier  was  still  in  great  evidence  and  it  was  still 
impossible  to  generate  a 
satisfactory  model  wave. 


Similarity 

Rules  for  Wave  Periods 


70.  In  all  three 
cases  studied  in  the 
flume  (undistorted  pro¬ 
file,  distorted  profile 
without  screens,  and  dis¬ 
torted  i  ofile  with 
screens),  the  resonant 
characteristics  of  the 
basin  interfered  with  the 
functioning  of  the  model. 
Recognizing  the  fact  that 
the  resonance  will  al¬ 
ways  be  strong  unless  ar. 
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Fig*  39*  Typical  wave  record  in 
distorted  flume  with  screens 


TERMINOLOGY:  £-1/  L 

i  =  FORESLOPE  LENGTH 
L  =  MODEL  LENGTH 
H  =  DEPTH  AT  WAVE  GENERATOR 
T  =  RESONANT  WAVE  PERIOD 
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Fig.  40.  A  comparison  of  theoretical  and  observed  periods 
for  different  model  configurations 

extremely  large  model  is  built,  it  was  realized  that  the  model  should  be 
constructed  to  resonate  at  the  desired  period. 

71.  A  theoretical  analysis  of  a  resonant  model  is  included  in 
Appendix  E.  The  results  of  the  theoretical  development  were  verified  by 
the  preliminary  model  tests  with  the  correlation  between  the  observed  and 
the  calculated  period  shewn  in  fig.  40.  The  theoretical  analysis  also  in¬ 
dicated  the  similarity  criteria  to  be  used  in  designing  a  resonant  model. 
The  ratio  of  model  period  to  prototype  period  was  shewn  to  be 
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(36) 
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rS 


1 

1/2 


where 


H 

r  =  distortion  ratio  =  — 

faV 

S  =  Sv  =  vertical  scale 
T  =  period  of  oscillation 


Furthermore,  the  resonant  period  for  a  triangular,  two-dimensional  basin 
with  a  closed  end  is  found 
to  be 

T  =  0.55  L 

I. 

o 

or  ,T  .  1/2 

(37) 


m 


T  =  0.55  .  „ 

m  Wi 


where 

L  =  length  of  the  basin 
a  -  bottom  slope 

72.  Employing  these 
criteria,  a  resonant  two- 
dimensional  model  can  be  de¬ 
signed  for  a  particular  ob¬ 
served  prototype  period.  A 
table  of  possible  model  di¬ 
mensions  is  presented  in 
fig.  4l.  These  values  were 
calculated  for  a  prototype 
period  of  28  min  and  a  bot¬ 
tom  slope  of  0.0060,  both 
of  which  conditions  were  realized 
tsunami. 


r  S 

Tm  =  0.55 
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r  =  DISTORTION  RATIO 


_  L0p  _ 
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Fig.  4l.  Model  lengths  for  given 
scales  and  distortions 

at  Crescent  City  during  the  1964 
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WAVE 

GENERATOR 


Problems  for  Future  Sti 


73.  One  of  the  main  findings  of  the  aforementioned  tests  is  that,  in 
the  model,  the  length  of  the  sea  area  between  the  wave  generators  and  the 
harbor  should  be  great  enough  to  create  at  least  one  cycle  of  the  wave  ap¬ 
proximating  that  of  the  prototype.  In  the  interest  of  economy  of  construc¬ 
tion  and  to  avoid  an  excessively  long  model,  it  is  necessary  that  the  model 
be  distorted.  A  further  requirement,  as  can  be  deduced  from  the  analysis 
of  scale,  is  that  the  model  be  constructed  affinely  with  the  prototype. 

This  would  mean  that  the  bottom  contours  of  the  area  from  the  harbor  to 
the  generator  should  be  sloping.  Based  on  these  requirements,  the  model 
dimensions  necessary  for  different  depth  scales  and  different  distortions 
are  summarized  in  fig.  4l.  A  slightly  modified  wave  generator  (shown  in 
fig.  42)  will  be  used  in  any  subsequent  testing. 

74.  Among  the  problems  to  be  considered  is  the  manner  of  controlling 
the  air  valves  that  would  lead  to  the  formation  of  a  well-shaped  wave.  The 
efficacy  of  screens  to  eliminate  the  secondary  undulations  related  to  the 

efflux  of  water  from 


J-»i.«NCXO  VALVM 


counter  to 
TO  VACUUM  MANIFOLD 


ftCMOVAaiC  THROAT  AUCTION 


the  generator  needs 
to  be  examined  anew. 
Another  quite  in¬ 
triguing  problem  is 
the  feasibility  of 
utilizing  the  waves 
repeated  in  the  flume 
following  the  initial 
operation  of  the  gen¬ 
erator  tank. 

75.  One  basic 
problem  remaining  is 
the  effect  of  distor¬ 
tion  on  the  runup. 

Ihe  necessity  of 
the  study  was 


Fig.  42.  Modified  tsunami  generator 


particularly  emphasized  by  Professor  R.  L.  Wiegel  in  a  conference  held  at 
the  WES  August  8-9,  1968.  In  the  development  of  the  analysis  for  the 
scales  it  was  shown,  tacitly,  thet  similarity  of  the  wave  quantities,  par¬ 
ticle  velocity  and  wave  elevation,  exists  in  the  flume  independently  of 
distortion,  provided  that  the  period  of  the  model  is  correctly  selected. 
This  is  a  consequence  of  the  dynamic  and  kinematic  equations  applicable  to 
the  flume.  The  argument  cannot  be  applied  to  the  runup.  Although  resort 
may  be  made  to  another  analysis  to  establish  the  runup  extensions,  it  would 
be  more  satisfactory  to  resolve  the  question  experimentally. 
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BAR?  VI:  C0BC1CSI08E 


76.  Based  on  the  results  of  the  material  presented  in  this  report) 
it  is  concluded  that: 

a.  The  nave  characteristics  that  should  be  used  in  the  proposed 
model  testing  program  to  simulate  the  highest  vac re  of  the 
I96I  tsunami  at  Crescent  City  are  a  23“*t,  sB-mln-period, 
rising  tide  with  a  motHm  runup  to  about  the  +20  ft  mllw 
contour. 

b.  Shelf  resonance  in  the  offshore  Crescent  City  area  is  a  key 
factor  in  the  unusual  height  and  duration  of  the  1964  tsu¬ 
nami,  and  probably,  but  to  a  lesser  extent,  of  all  other 
tsunamis  causing  a  water- surface  rise  at  Crescent  City. 

c.  Linear  spectral  analyses  performed  on  a  portion  of  the  i960 
tsunami  record  at  Crescent  City  show  that  the  contained 
periods  of  this  tsunami  coincide  with  those  derived  by 
Keulegan  and  Wilson  to  be  the  resonant  oscillations  of  the 
offshore  shelf  area. 

q.  A  digital  computer  program  that  allowed  for  the  distortion 
present  in  maps  drewn  by  Mercator  projection  was  coded  for 
tsunard  refraction  across  the  Pacific  Ocean  area  and  was 
tested  successfully.  Refraction  diagrams  computed  and 
drawn  by  the  above  program  for  three  recent  tsunamis 
(Alaska,  1964;  Chile,  i960;  and  Kamchatka,  1952)  estab¬ 
lished  nearshore  wave- front  patterns  at  Crescent  City. 

These  nearshore  wave  fronts  will  dictate  the  wave  gener¬ 
ator  alignment  in  any  subsequent  model  testing  program. 

The  refraction  diagrams  also  showed  that  no  overall  tsunami 
height  increase  was  caused  by  refraction  in  the  open  ocean. 

e.  The  probability  of  occurrence  studies  indicated  that  tsunami 
heights  at  Hilo  and  Crescent  City  are  similar  and  that  the 
1964  tsunami  was  a  singular  event  with  only  a  20  percent 
probability  of  occurrence  during  any  45-year  period. 

f .  The  use  cf  an  undistorted  tsunami  model  of  the  Crescent  City 
area  will  not  be  feasible;  and  the  sea  environment  of  the 
model  should  be  of  such  a  magnitude  that  the  fundamental 
period  of  the  resonant  oscillations  of  the  model  sea-to- 
land  area  corresponds  to  the  largest  period  of  the  selected 
test  wave. 

Additional  wave-flume  tests  are  required  to  determine  the 
optimum  linear  scales  of  the  model  and  the  required  sea  area 
that  must  be  reproduced  to  ensure  accurate  reproduction  of 
prototype  tsunami  action  in  the  harbor  and  overbank  areas  at 
Crescent  City. 
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Table  3 


Comparison  of  Honolulu  and  Hilo  Heights  for  Tsunamis 
with  Large  Runup  at  Hilo 


No. 

Date 

Generating 

Area 

Hilo,  ft 

Hono¬ 
lulu,  ft 

Hilo 

Reduced,*  ft 

Crescent 
City,  ft 

2 

1896  June  15 

Sanriku 

7-9 

0.33 

0.90 

1.81 

3 

1906  Aug  17 

Talcahuano 

4.9 

0.33 

0.90 

0.89 

4 

1922  Nov  10 

N.  Chile 

6.9 

O.98 

2.68 

1.84 

5 

1923  Feb  4 

E.  Kam¬ 
chatka 

20.0 

2.9 

7.91 

0.9 

* 


Values  shown  for  "Hilo  Reduced"  are  Honolulu  values  multiplied  by  2.73- 


APPENDIX  A:  COMPUTER  PROGRAM  FOR  TSUNAMI  REFRACTION 


Lanatlon  of  Variables  and  Program  Listing 


1.  Tbe  program,  for  which  a  listing  is  given,  is  written  in  FORTRAN 
II  and  has  been  run  successfully  on  the  GE-225  computer  at  the  U.  S.  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  Miss.  The  program  gives 
output  suitable  for  the  off-line  Calcomp  750  plotter  system  at  the  same 
installation. 

2.  The  compilation  takes  approximately  150  sec,  and  the  program 
executes  at  the  approximate  rate  of  three  points  per  second  on  the  GE-225. 
These  execution  times  include  reading  the  data  for  a  35  x  35  grid  and 
printing  output  every  50  points. 


3.  Input  Variables 


I'M . The  number  of  grid  intersections  along  the  x-axis. 

NN  .  .  .  . . The  number  of  grid  intersections  along  the  y-axis . 

NOR . .  The  number  of  rays  to  be  refracted  during  one  run. 

NFRINT  .  Frequency  for  printed  output. 

GRID . The  length  in  inches  along  one  edge  of  a  grid 

square . 

H(l,J) . Depth  data  at  the  grid  intersections  (Note: 

I  as  MM,  J  £  NN) . 

X,Y  .  Coordinates  (in  inches)  of  the  starting  point  for 

j  each  ray. 

ANGIE . Initial  direction  of  the  wave  ray,  measured  in 

degrees  counterclockwise  from  the  positive  x-axis . 

DELTAX . Computational  x  increment  (&x)  measured  in  feet. 

BIST . The  number  of  feet  corresponding  to  the  length  of 

a  grid  square  (measured  at  the  equator). 

CCNV . Conversion  factor  to  change  the  depth  readings  to 

feet . 

TIMEX . Dictates  the  number  of  seconds  of  actual  wave 

travel  between  tic  marks  placed  on  the  ray. 

(Note:  The  wave-front  positions  may  be  marked 
TIMEX  seconds  apart  on  the  ray.) 


TIMS . The  tine  in  seconds  at  the  beginning  of  a  wave 

ray  confutation. 

U.  Output  Variables 


!  As  previously  described  except  allied  to  the 

X,Y  >  .  . . successive  computational  steps  instead  of  the 

»  initial  conditions. 

DEPTH . The  interpolated  depth  in  feet  at  successive 

computational  points. 


5-  Variables  hot  Previously  Defined 

HE( l) ,HE(  2) , . . ,HE( 10) .  .  Depths  at  grid  intersections  used  for  surface 

fitting  computations. 

SMALL* ;SMALLY  .  Distances,  in  feet,  of  the  point  of  interest 

measured  fro.*  H(I,J)  . 

3H  3H 

PHWRTXjPHWRTY  ....  found  from  the  surface  fitting  procedure . 


D2YDX2 . £4 

dx 

am  ....  . . & 

DELTA’/  .........  The  increment  of  y,(Ay)  ,  computed  from  a  given 

increment  of  x,(Ax)  . 

DELTAT . The  tinu  Increment  between  successive  points  of 

calculations . 


6.  It  must  be  noted  that  the  format  for  the  depth  data  is  not  vari¬ 
able  but,  of  course,  may  be  changed  to  suit  the  needed  requirements. 

7.  The  plotting  steps  are  incorporated  as  an  integral  part  of  the 
program.  A  set  of  x-,  y-axes  are  drawn  with  MM  tic  marks  placed  on  the 
x-axis  and  HN  tic  marks  placed  on  the  y-axis.  These  tic  marks  are  placed 
GRID  inches  apart  and  the  axes  form  a  reference  by  which  the  plotter  output 
may  be  overlaid  on  the  original  map  and  the  rays  (which  are  drawn  by  the 
Calconp  plotter)  traced  from  the  plotter  sheet  to  the  original  map. 

8.  Because  of  main  memory  limitations  in  the  GE-225,  the  maximum 
aUe  grid  that  can  be  used  is  35  x  35u 

9.  A  progrea  listing  and  an  example  of  input  for  the  computer  pro¬ 
gram  are  found  on  the  following  pages. 


PBOGBAM  LISTING 


_ c  THEORETICAL  TSUNAMI  STUDY  FOR  CRE8CENT  CITY. CALIFORNIA _ 

C  JOHN  HARRISON  2I-Z5-074 

C  PROGRAM  computes  and  plots  refraction  diagrams  for  a  tsunami 

DIMENSION  H(35«35)«HE(101 

TANF I AN6RAD 1  »S1  NF I ANSRAPJ /COST  I  AjjGRAD ) _ 

0«32. 174 

READ  1,HN,NN,NQR.NPRINT.6RID _ 

1  FORMAT  (4H0.F10 .5) 

_ LL»0 _ 

READ  2*(IM(I»J| , I*1,MM| , jal.NNI 

_ _ 2_FQRHATJ.7fiP.ttJ _ 

CALL  PLOTIO. 0,-30. 0,-31 

_ —  CALL  fLflIU).0tS,.0i=3l _ 

MMP«  MM-1 


_ MMPJ-hMrJ _ 

call  tmkio.o»o.o,8Rio»mmp< 0.1,0.0.1.0*1.01 

_ CALL  TMKI0.0»0.0,6B1P,NNP. -0.1, <0.0, 1.0, 1.01  _ 

CALL  PLOTIO. 0.0. 0,-31 

_ -1QQ2  READ  3tlf.lt AMftLfjJfiLIAJUQiAlxCQMYi  T 1  H5X/.T.IHE _ 

MPRINT  ■  1 

_ -  ..mu.fl _ 

3  FORMAT ISFIO .31 

_ HEAll _ : _ 

tihes«timex 

_ PRJHT-JL4 _ 

14  FORMAT (1H1«2X«10H£NGLE(OE61 .5X.16HX  COORDINATE; IN) ,5X,16HY  COORDIN 

_ lATEHNI.OK.OHTlMEISEC;. 12X.10HDEPTmrr.il _ 

999  ANGRAD«ANBlE*0. 0174532929 

.  J02.JL»JLt±l(/6m _ 

J»1.*Y/GRID 

. -025.HE  Ii»  «HLL»Jl»CQHy _ - _ _ 

HE (21 "Ml  1*1, J!*CONV 

_ HE131»H(1*2,  Jl«-CONV _ 

HE(41«H( I»J*ll*CONV 

_ HE  L5J  _ 

HE  1 61 "Hi 1*1, J*1 J*CONV 

.  _ He  .171 1HJLI  Ilf  J-lL»CaNj> _ 

HE  1 81  «H  I  I»J-U*CONV 

- . . HE(91«Hn*g,J*ll«CQNV _ 

HE(10l«M  Wl.J»21*KONV 

. . 5.9.2 -XJLlI:^ _ 

yy»j-i 


XFEETiJBI  SJ_*X/MJ.tL_ 
YFEET"  DIST*Y/GRID 


SMALLX  «  <FEET-XX*D1ST 
SMALLY  •  XFEST-YY*DIST 


A3 


LISTING 


Cl* (HE 161-HE I 7!-2.*Bl*DJSTj/ (2. *0181*018X1 


nLi^nTBlijUDi!  imtitwj)  iioh  it  I 

jrgnmi 

graGSHTlM 

rcrnfjM 

™»i'!!5.}  3r»-| 

mmm 

^TnlMT7.H'a  ujftAO-rcf  'o  itt.ti  hi 


•MEUI-D2* 


C2*IH£ t91-HEill’32.*A2*D!ST-B2*DIST»4.*D2*D18T#DIST-E2*0lST*Dl8Ti7 
.*DIS7*DIST1 


A*|Al*A2l/2. 


C*CCl*C2l/2. 


E*|El*E2)/2. 


pmwrt ybB*c*sm all x*2 • *e*sm ally 


ni*ii«nn  nv)j 


60  TO  1001 


BJiyiTxo7i!iTn!^!crTir^7TrCl74T^g  a*gan 

r  V  T!T 1 1 V  1  rT  1*1  yrWMBBBB 


4004  MPRINT*! 

400iL PRINT  890*AH6LE.X.Y.TlHE, DEPTH _ 

000  F0RMAY;3X»F1 U«  4.9X»F7.3»l3X»F? .3»8X*F14.4»10X,F14 .31 

400;  60  TO  1814.8151. NEB _ 

614  CALL  PlOTiX.Y.3l 


AN8N6 W* AT ARF <  DYDX 1*97. 20977991 


DE|.TAT*(S«RTFIDELTAX*0ELTAX*0ELTAY*DELTAYI)/SQRTFI6*DKPTH| 
.4 


4001  HPR!NT*HPR!NT*1 


60  TO  816 


816  lFITIME-O. 01817. 818. 817 
ft 

817  IFIT1HES-TIMEI887.819.825 


1 1  wrin  i.ui  tJifi'frnn!!*  xift.ilki  mnn  ni  i 


SO  TO  888 

J387  TID1FF*TIHE-TIHE3 _ _ _ _ 

PL0TX*|S0RTF I6*DEPTHJ *T1D*FF *6RID1/ IS0RTF 11. *DYDX*DYDX)*DIST1 


m  WffiSBMm 


plcty*y 


PROGRAM  LISTING  (flont'd) 


CALL  TMKl»L0TX,PL0TY, 0.0, 1,0.1.  ANGLE.  1.0,1.01 


'TT.MI'1=!Mt 

■awflTgarri 

DTrSBtSiSjifll 

In  „  . 


oeltay«deltay*grio/dist 

JSMiBSUlAI _ 

y*y*delta» 


1.0,1.51 


.u  *:ou-i.:r 

irifi.T-in  rw 


inANGLE*86. 01 3002,3002, 3001 


3002  PRINT  3003 

3003  FORMAT I52HAN6LE  IS  BETWEEN  69  DEG  AND  271  DEG.  START  A  NEW  RAY1 
80  TO  1001 


rm'M’infj.a  'imio  'i  i  lfxi.i 

nn 

im 

wnnmrnm.inp. 

i  i7rrrrgrrTTTT7rTTn  > -  i  * . 

PTfMyfrfnrToiWBPBTHaTiaarirTTnPFiTBBTiT 


xxx«hm 

YYY«NN 


XXXX*(XXX*2. I*6RID 


1099  PRINT  5000 


60  TO  1001 

1000  iriYYYY-Y»S001.50l 
$001  PRINT  5002 


JJ I  rl-l  iJB  L'lm'imd 

fTrumTJfTVrai  U't/niiMiii 


$003  PRINT  5004 

-J^H0.4-flUSfeAT  1 31MRAY  IS  QQJ.NtJM 

1001  If;LL-NOR*102S. 1003. 1003 


PRINT  600  » ANGLE. X.Y, TINE. DEPTH 

_ *ElUNJL-3 _ 

READ  1,  JOHN 

. . -JNDL _ 


jggjffiLiwaL 


4*98. 

•852. 

4>*r. . 

ii%57 

4353. 

il7<n 

— 

3944. 

412*.  _ 

4418. 

AM*.  _ 

4*99. 

45*5. 

4571. 

4m. 

4*3*. 

4*>4. 

4?fS. 

4*9*. 

49*9. 

4**4. 

457*. 

4**5. 

4882. 

4123. 

4525. 

4553. 

4923. 

4249. 

4*Jl. 

47**” 

4**2. 

47*2. 

476*. 

4>«*. 

492*. 

4*M. 

45*1. 

4*93. 

45*7. 

4522. 

4233. 

4218. 

«29». 

44*8. 

4712. 

4753. 

4898. 

4751. 

422*. 

499*. 

5178. 

592*. 

595*. 

5168. 

4*11. 

45*2. 

3943. 

3*72. 

4342. 

453*. 

5998. 

5*45. 

4933. 

493J. 

4933. 

4933. 

4933. 

«933. 

4*33. 

5132. 

5*49. 

55IT 

KT. 

5443. 

*3*4. 

53*4. 

49*9. 

4992. 

5817. 

5344. 

54*2. 

551*  , 

5558. 

5*84. 

5*51. 

55*9. 

5*97. 

Ml*. 

5946. 

417*. 

445*. 

4582. 

4*43j_ 

4129, 

5126. 

57*3. 

5857. 

5*29. 

5U2. 

3118. 

5116. 

511*. 

5118. 

5221. 

5233. 

5192. 

3144. 

5*41. 

4*7», 

4787. 

4922. 

4434. 

4694. 

4779. 

4754. 

«8^. 

4*58. 

4558. 

4533. 

4582, 

S321. 

5379. 

4751. 

4523. 

-  45*2. 

4774. 

4696.' 

4*1*. 

513*. 

5*2*. 

58**. 

5241. 

523*. 

5293. 

5297. 

5297. 

5379. 

5362. 

e 5235. 

543*. 

5399. 

3351. 

5298. 

5135. 

521*. 

8  **94. 

4957. 

48*9. 

4944. 

5928. 

517*. 

5192. 

1  511*. 

5291. 

493*. 

4*97. 

5153. 

47*5. 

4496. 

2  45*3. 

5153. 

5171. 

4971. 

5398. 

521*. 

5281. 

5  51*9. 

599?. 

5414. 

9441. 

538*. 

549*. 

5369. 

sc  5*15. 

5599. 

5578. 

554*. 

54*9. 

5347. 

5349. 

|  5298. 

5922. 

4928. 

3982. 

51*2. 

5313. 

5228. 

S  52*6. 

5248. 

4923. 

4942. 

5112. 

499*. 

4967. 

5917. 

5518. 

55*8. 

5663. 

5*92. 

5574. 

5567. 

5459. 

538*. 

5*95, 

55*1. 

5483. 

5591. 

5**5. 

5*95. 

5581. 

5*42. 

5598. 

5477, 

$349. 

5591. 

5265^ 

5553, 

493*. 

4845. 

5124. 

5186. 

5248. 

5378. 

5311. 

5234. 

525*. 

5377. 

5361. 

5242. 

5391. 

5655. 

6994, 

4251.  _ 

5890. 

5765. 

57«7. 

5935. 

5990. 

5745. 

5725. 

56<*. 

5646. 

5646. 

5*2*  ... 

562*. 

5622. 

5*17. 

5546. 

5549. 

5651. 

5112. 

5*73. 

51*8. 

5112. 

5252. 

$269. 

5126. 

522*jl  _ 

_  5lg8. 

_ _ 5255. 

5555. 

5521. 

559*. 

5603. 

5793. 

6182. 

6224. 

6192, 

6165. 

6995. 

6049. 

5973. 

5795. 

5598.  _ 

5579,  _ 

5747. 

5717. 

5651. 

5722. 

5644. 

5619. 

5642. 

5599. 

5540. 

5695. 

5514,  . 

54*2, 

5373. 

5952. 

4787. 

4897. 

4862. 

5253, 

5231. 

5362. 

5593. 

5630. 

5328. 

5766. 

5903. 

6171. 

6250. 

6129. 

5960. 

6239. 

6103. 

6991. 

593*. 

58*2. 

5678. 

5904. 

6154, 

6464. 

5544. 

5544. 

5649. 

5649. 

5649. 

5649. 

5544. 

5*96. 

5620. 

5630. 

5606. 

5681. 

5361. 

5*43. 

5435. 

5415. 

4*52. 

4970. 

5504. 

5561. 

5844. 

5925. 

6278. 

6252. 

6178. 

6242. 

6257. 

6112. 

5968. 

583*. 

6011. 

6208. 

6256. 

6347. 

6332. 

pwrwKoy’O 


ABBBMUX  ft  HMKFOMttlO*  Of  WV8  WTSM  01  IB  BARM’S  SURFACE 
®  A  MDBCA3DR  EEFHB99RAH0M  OF  XRAT  SURFACE 


1.  Latitudinal  and  meridional  lines  of  &  sphere  nay  be  represented 
an  %  mp  so  that  tine  lstitadM  are  horizontal,  parallel  lines  and  the 
longitudes  are  vertical,  parallel  lines.  Such  a  representation  is  tenaed 
a  Mercator  projection  and  is  arthoaorphic  (i.e. ,  the  distortion  is  unifora 
in  all  directions  at  each  point  on  the  aap). 


2.  Let  Axq  *  eleaentary  length  on  a  parallel  of  latitude,  t 

dgrQ  *  eleaentary  length  on  a  aeridian,  0 

R  ■  radius  of  the  sphere 

Ax  *  cos  f 
o 

Ayo«  na# 

3.  Let  Ax  and  Ay  he  corresponding  lengths  on  a  Mercator  aap. 
Since  the  projection  is  orthaaorphic. 


Ax  Ay 

Since  two  consecutive  Meridians  separated  by  A^  are  represented  on  the 
Mercator  aap  by  parallel  lines  separated  by  Ax  »  R ,  we  have 


WaSTi  *  5ec  ♦"" 


Therefore, 


«  n  »  sec 


»  R  sec  f 


y  *  R 


\J  sec  f  df 


y  «  R  log  tan  ( J  +  |) 


B1 


/ 

% 

Considering  a  unit  sphere  where  R  =  1  , 

er=tan(j  +  i)  (Bl) 

Also,  it  can  he  shown  that 

ey  =  sec  Jr  +  tan  Jr  (B2) 

4.  The  above  derivation  shows  that  the  distances  on  a  sphere  are 
magnified  by  the  factor  sec  ijr  on  a  Mercator  representation  of  this  sphere. 
The  present  problem  is  to  determine  a  method  of  plotting  wave  motion  which 
takes  place  on  the  earth's  surface  on  a  Mercator  map.  This  is  accomplished 
in  the  present  methodology  by  multiplying  the  actual  input  depths  by  the 
secant  squared  of  the  latitudinal  angle  Jr  .  Since 


and 


V  =  ygHQ  sec2  Jr  =  sec  ♦  ^gH~ 

V  =  V  sec  Jr 

o 


The  time  element  dt  is  invariable  since 


dt 


ds 

V 


dS  sec  Jr 
o _ _ 

V  sec  * 
o  T 


5.  This  subjective  reasoning  can  be  proved  in  more  definite  form. 

Referring  to  fig.  Bl,  the  basic  refraction  equation  is: 

✓ 


d0 

ds 


sin  9 


ac 

- 


cos  8 


ac  n 
^ ) 


where 


sin  6 


ds 


_ 

V^(dx)2  +  (dy)2 


B2 


Now 

ey  »  sec  ♦  +  tan  t 


and 

so 


tan  f  a  ey  -  sec  f 
tan2  f  a  e2y  -  2ey  sec  f  +  sec2  t 
sec2  t  -  1  °  tan2  ♦ 

2ey  sec  f  »  X  +  e2^ 
sec  ♦  “  |  (ey  -  e“y)  -  cosh  y 


The  depths  are  transformed  so  that 

H  a  Hq  sec2  f 


(B4) 


In  a  sea  of  uniform  depth 

H  ■  H0  sec2  ♦  ■  Hq  cosh2  y 
B4 


tanh  y  dy 


Integrating 


or 


rV-2 

T)  +  Tf 

I51  -  *  2  tanh  y  dy 

in  T)  -  in  (l  +  Ti)  =  2  in  cosh  y  +  in 


7)  *2 

rhi  ’  ci cosi  y 


(B9) 


since 


\  -  (I)'  ■ cot2  ♦ 


where  f  is  the  angle  between  the  ray  and  the  x-axis.  At  y  =  0  ,  is 
evaluated  as 

.2  \ 


1  1  + 


6.  For  exanqale,  if  the  ray  follows  the  equator,  0=0,  cot  0  =  »  , 

and  C,  becomes  unity.  If  the  ray  follows  a  meridian.  0  =  90° » 

2  1  2 
cot  0  »  0  ,  and  ^  »  0  .  Now 

T)  =  C2  cosh2  y  +  T|C2  cosh2  y 


or 


U 


2  2 
C1  cosh  y 


dx 

3y 


~ 2 - 5“ 

1  - 

cosh  y 

ta\2 

cosh2  y 

dy/ 

“  T2 — - 2 — 

-  cosh  y 

cosh  y 

1  -  sinh  y 


B 6 


it'tn. 


laf  «*******■#> 


Then 


Combining  equations  BI3  and  BlH  into  equation  BI2  one  obtains 


&g  _  cosh  y 
**  Vk2  -  sin2  y 


(B15) 


Tho  tstnad  travel  tine  can  be  represented  by  the  equation 


Where 


is  the  path  of  the  ray 

V  is  the  velocity  over  the  path  on  a  Mercator  nap 


The  velocities  have  been  assuned  to  be 


(B16) 


Therefore 


v0  sec  ♦ 

.  1  r  ds 

VoJ, 


(B17) 


7.  In  order  to  test  the  depth  transformation ,  consider  a  vave  ini¬ 
tiated  at  x  -  0  ,  y  =  0  ,  and  traveling  to  the  aeridian  x  =  ^  .  Equation 
B17  becanes 

,  r"A  \h  a* 


ns; 

cosh  y 


Equation  B15  becenes 


\2  K2  -  sirih2 


cosh  y 


K2  -  sinh2  y  +  cosh2  y  1  +  K2 


cosh  y 


cosh  y 


(Bl6) 


Therefore 


t  = 


/tt/2 

dx 


- g— 

o  cosh  y 


(B19) 


Proa  equation  Bll,  equation  B19  becomes 

-tt/2 


t  = 


and,  when  xq  =  0  , 


Vl  +  K2 


/TT/C 

1  + 


dx 


„2  .  2  ,  N 
K  sin  (x  -  xq) 


t  a 


I -  fn/2 

Vl  +  K2  |  dx 


•J  IT 


K2  sin2  x 


(B20) 


Equation  B20  can  be  integrated  by  the  formula 
dx  1 


/?T 


”5 - 2 - 7== 

b  sin  x  •/  2  ,2 

aVa  +  b 


tan 


If  Y  a2  +  b2  tan  x 


so  that 


t  = 


Vl  +  K2 


. -  tan"1  (  Vl  +  K2  tan  x  ) 


LVi  +  k2 

1  /IT  TT 

*  "  V  (2  "  °)  “  2V 


tt/2 


8.  It  can  be  concluded  from  the  development  given  above  that  all 
waves  starting  from  (x,y)  =0,0  reach  the  meridian  x  =  ~  at  the  same 
time.  Since  this  fact  can  be  immediately  understood  to  occur  on  a  sphere 
covered  with  a  uniform  depth  of  ocean,  and  has  been  shown  to  occur  under 
the  proposed  transformation  of  depths  onto  a  Mercator  map,  the  proposed 
transformation  is  justified. 

9.  A  refraction  diagram  drawn  by  digital  computer  techniques  uti¬ 
lizing  the  above  presentation  (and  previously  discussed  at  length)  is  shown 
in  fig.  B2.  This  refraction  takes  place  across  an  idealized  ocean  which 
everywhere  covers  the  earth’s  surface  to  an  equal  depth. 


B9 


DEGREES  LONGITUDE 


Fig-  B2.  Representation  of  wave  refraction  on  Mercator  map 


A2J3NDIX  C:  TSUNAMI  REFRACTION  IN  THE  AREA 
OF  A  1HEOHETICAL  ISLAND 


and  free  equations  C5  and  C6 

_ r  +  2r'2  -  rr" _ 

(p2  ♦  r,2)(r'  cos  6  -  r  sin  6) 

Tbe  tsne  of  equation  C7  are  expanded  as  follows: 


a*  ar  ax  a®  a* 

*5,  figai  +  Mas 
3j  dr  ay  ad  ay 

and 


(C8) 


=  cos  8  ,  ^  =  sin  6 

3*  ay 

jg  _  sin  6  38  _  cos  8 

a*  =  *  r  ay  r 

Bow  since 

dy  _  r*  sin  8  ^  r  cos  8 
dx  r*  cos  8  -  r  sin  6 


2H \dx  ax  ay/ 


j  (r*  sin  8  +  r  cos  6)  -  (r*  cos  0  -  r  sin  8) 

_ I _ _ _ 0* _ I _ ayj 


r '  cos  9  -  r  sin  6 


ayJ  (c?) 


Considering  tnese  terns  separately 


(r*  sir.  0  +  r  cos  e)  ~  =  |r'  sin  0  +  r  cos  dj  ^  c.os  0  ^  -  — y  —  j 


=  r'  sin  8  cos  0  ^  +  r  cos2  8  ~  —  sin2  0  -  sin  0  cos  0  ^ 

ar  r  06  3“ 


C3 


t-  -  KXlfl.rir^  ♦*»* 


(r*  cos  0  -  r  sin  0)  =  (r*  cos  0  -  r  sin  0 j  ^ sin  0  +  cos  0  j 

=  r ’  cjs  0  sin  0  -  r  sin2  0  ^  +  ~  cos2  0  ^  -  cos  0  sin  0 

dr  oT  r  00  00 

Therefore 

(r’  sin  0  +  r  cos  0)  ^  -  (r*  cos  0  -  r  sin  d)  ^  =  r  ^  ^  ^ 

0*  oJ  or  r  o9 

and 

i  (ras _£i$n 

1_  (  &  35  .  M\  =  gH  V  ar  r  ae  / 

2H  V  dx  3x  fcy/  r*  cos  0  -  r  sin  0 

Now,  combining  equations  C7j  C8,  and  C9  yields 


o  2 
r~  +  2r* 


Equation  CIO  is  the  refraction  equation  in  polar  coordinates. 
3.  Let  us  assume  that  the  island  is  symmetrical  so  that 


(CIO) 


3&  = 
m 


0 


Therefore 


2 

r 


+  2r 


,2 


(cn) 


Tne  hypothetical  island  is  illustrated  in  fig.  Cl  and  is  defined  by 
the  equations 


C4 


a  <  r  <  b 


The  solution  of  equation  Cl4  is 


8  =  Cx  log  r  i-  C2  (C15) 

In  order  to  solve  for  the  two  constants  in  equation  C15,  consider  a  plan 
view  of  the  given  island  shape  (fig.  C2).  The  x-axis  is  situated  parallel 


y 


Fig,  C2.  Flan  view  of  theoretical,  island  showing 
refraction  of  one  wave  ray 

to  the  incident  ray  direction  and  all  angles  are  measured  in  the  usual 
co'.uiterclocjkwise  direction  from  the  positive  x-axis. 

4.  At  the  point  of  entrance  of  the  ray  to  the  island's  influence 

Mrt“K-i)*-cot,b  (cl6) 

Therefore 


dr  r. 

o 


C6 


From  equation  C15 


8„  =  135°  =  3? 


tan  6^  =  -1 

so  that 

6a  -  $ +  log  t 

Typical  values  of  0&  (with  6^  =  for  various  values  of  the  ratio 
are  tabulated  below: 


a 

b 

log  | 

0a 

0.5 

-0.693 

1.666 

0.4 

-0.916 

1.443 

0.3 

-1.204 

1.155 

0.2 

-1.609 

0.750 

0.1 

-2.303 

0.057 

6.  Fig.  C3  slows  a  typical  refraction  pattern  around  the  given  island 
shape  to  which  the  specific  dimensions  shown  in  the  figure  have  been 
assigned.  The  results  given  by  the  accepted  computer  program  agree  very 
closely  with  those  produced  by  the  theoretical  treatment. 


a'lfB 


Fig.  C3.  Typical  refraction  pattern  around  the 


APPENDIX  D:  SHELF  OSCILLATIONS  IN  THE  CRESCENT  CITY  AREA 


1.  One  may  consider  the  oscillations  of  an  ocean  shelf  by  analyzing 
the  simplified  two-dimensional  model  shown  in  fig.  Dl.  The  model  config¬ 
uration  consists  of 


Fig.  Dl.  Two-dimensional  model  shelf 


a  uniformly  sloping  shelf  which  drops  off  to  infinity  at  the  ocean  end  and 
is  terminated  by  a  vertical  wall  near  the  land. 

2.  Let  the  distance  x  be  measured  from  the  point  of  intersection 
of  the  still-water  surface  and  the  shelf  slope.  The  shelf  then  extends 
from  the  point  x  =  x^  to  x  =  L  .  The  depth  of  water  varies  uniformly 
as 


x 

o  L 


(Dl) 


3.  Denoting  the  surface  disturbance  as  h  ,  the  velocity  of  water 
in  a  section  x  by  u  ,  and  the  time  by  t  ,  the  equation  of  motion  be¬ 
comes 

Bu  3h  (D2) 

at  “s  ax 

where  g  is  the  gravitational  acceleration.  The  conservation  of  mass 
equation  is 


Preceding  page  blank 


5.  A  second  relation  similar  to  equation  D17  inroting,  llk**i*ar 
tLe  constants  A  and  B  should  be  obtained.  At  the  jriut  f  *  §, 
particle  velocity  is  defined. 

t.  te-^xssider  equation  332  avfl  KjfaetiCute  from  equation*  D5  end  C6. 


c'u  *  -£ 


*5 


(318) 


ou 

L  45 


i©i9> 


This  cas  also  be  uritles  as 


- 


(5») 


5ifftredUtii^  equation  Di5  with  reipeci  tp  f  , 


ard  sicce 


and 


d*  .  •  : .  , 

—  -  Aj  .xi'_  4.  SK  ,9, 

28  o  o 


j;(e)  =  -Jx(e) 


k;(6)  =  -1^(0) 


we  nave 


|3  =  -  AJ^e)  -  BK^Ce) 


(D21) 


Now 


dH  .  dji  d9  and  d0  _  2K 
d?  ~  d9  d5  anQ  d?  ~  0 


D4 


sc  that 


! 


Y 

*3  a  « 

S'e  d6 

(D22) 

L  ! 

i  i 

r  i 

i'  ■ 

Substituting  equation  D22  into  equation  ”£1  yields 

5?  -  f  ["!<*>  *  *!<•>] 

t  : 

i 

[ 

and  fit*  equation  D20  the  Telocity  relation  ii 

r 

\ 

\  ; 

L  I 

r  i 

r  i 

t  i 

*o  •  f  ["l(9>  * 

where 

(323' 

l  : 

*  J 

I  : 

b 

IT; 

.  =  ^  ^ 

k 

t 

i 

t 

At  t*str  coast,  f  *  5^  ,  thi  jfartiils  velocity  u  -  C*  ,  and  hence  V 
«ol 

vanishes 

■ 

AX.it.  1  -  SJ^O^-  -  0 

(set) 

t 

where 

0,  .25V2  ^ 

(D25) 

t 

Solving  for  A  and  B  from  equations  Dl8  and  D24 

V«i>\ 

A‘  ^ 

(326) 

> 

,  B  "  A 

(D27) 

| 

i, 

where 

‘ 

F 

A  -  Jo(WV  -  Jl(9i>Ko<eL> 

(D28) 

r 

1  Substitute  into  equation  D15  so  that 

D5 


Table 


*  Values  for  V&riou 


0 

*>.226 

2.276 

0.01 

5.12$ 

2.202 

0.062 

5.0.13 

2.112 

0.04 

4.880 

2.016 

0.09 

4.511 

1.797 

0.l6 

4.04$ 

1.35& 

Note:  N  = 


lift,  w 

L  2K1/2 


APPENDIX  E:  STANDING  WAVES  IN  TSUNAMI  MODEL  CHANNELS 


I.  For  simplicity  in  argument  it  will  be  assumed  that  the  harbor 
and  bay  environments  in  the  tsunami  model  can  be  approximated  by  the 
special  channel  configuration  shown  in  fig.  El. 


Along  the  distance  QA,  the  harbor  and  bay  area  is  correctly  contoured. 
Along  the  distance  AB,  the  water  depth  is  con*t:'"t.  The  wave  generator 
placed  at  B  constitutes  a  barrier,  and  at  this  point  the  particle  velocity 
u  vanishes  when  flow  from  or  into  the  generator  has  ceased.  For  the  suc¬ 
cessful  operation  of  the  model  it  is  desired  that  the  period  of  the  stand¬ 
ing  waves  in  this  particular  channel  be  greater  than  or  equal  to  the 
tsunami  model  period.  Although  in  principle  it  is  always  desirable  that 
harbor  and  bay  contours  be  reproduced  up  to  the  foot  of  the  generator, 
this  may  not  be  possible  because  of  cost  considerations.  In  any  case  one 
must  have  an  expression  for  the  period  of  the  slowest  oscillation. 

Law  of  depth: 

H  -  H  7  0  <  x  <  i 

o  i 

H  =  H  i  <  x  <  L 

o 

Dynamic  equation: 

&  _  _K  & 

St  "  s  &x 


(El) 


(E2) 


El 


Preceding  page  blank 


Continuity  condition: 


l*i<H“>  =  0 


Details  of  these  relations  can  be  found  in  reference  27.  From  equations 
E2  and  E3 


Introducing 


equation  E4  gives 


Introducing 


h  =  7]  cos  at  ,  cr  = 


a2T]  +  g|E(H|3)=0 


,  _  x  -  _  J 

5  "  L  »  51  ~  L 


the  last  equation  becomes 


whereupon 


One  must  consider  two  solutions  of  the  equation,  first  for  the  range 


0  <  5  ^  and  second  for  the  range  <  §  £  1  .  For  the  first  range 


H  §, 

O  31 


(E10) 


and  equation  E9  becomes 


d  T]  dT| 

5  +  dT +  Kq\ 


(Ell) 


Subject  to  the  condition  that  T]-,  is  finite  at  f  =  0  }  the  proper 

27  • 

solution  is  (Lamb  ) 


\  -  ai  jc  ( 251/2  5i/2  k1/2) 


(E12) 


For  the  range  §  <  |  £  1  ,  since  jj-  =  1  ,  equation  E9  becomes 

o 


•  — jr  +  K7]  =  0 

dt 


and  the  solution  is 


\  =  A2  cos  K1/2  §  +  B2  sin  Kly/2  § 


(E13) 


(El4) 


Now  at  §  =  1  ,  u  =  0  and  thus  ^|-  =  0  .  Ey  differentiation 


=  K1/2  (  -Ag  sin  K1/2  5  +  Bg  cos  K1/2  § ) 


(E15) 


A2  sin  W2  =  B2  cos  K1/2 


X 


01 


B2  =  A2  tan  K1/2 


(El6) 


At  the  section  through  A,  5  =  both  T)  and  ^  are  continuous. 

Hen  jo, 


and 


(E17) 


<%  dTl2 

W 

Then  from  equations  E12  and  El4 


d%  ~  dg  >  ?  =  5X 


J 


liJo(£Vl/s)= 


Ag  cos  +  B2  sin  ¥?'  c  ^ 


K1/2 


(E18) 


and 


Al^2  Jo  (  25!k1//2  )  =  (  -a2  sin  ^  +  B2  cos  K1/2  ^ 


or 


(■V0 


K1/2  C,  +  Brt  cos  K1/2 


-  -Ag  sin  IT'  -  5X  +  Bg  cos  I C'  ^ 


(E19) 


v; 

Making  use  of  equation  El6 

vJ^W2)  =  Ag  (cos  K1/2  5^  +  tan  K1/2  sin  1(1/2  0  (E20) 


and 


v;  (*/'■)•■,(  -  sin  K1/2  51  +  tan  K1/2  cos  K^^2  5-^j  (E21) 


E4 


l 


Dividing  equation  E21  by  equation  E20  and  utilizing  the  rule 


(  )  -  -Jx(  )  (E22) 

one  has 


The  roots  of  this  equation  of  least  value  for  specified  values  of  , 
through  equation  E8,  determine  the  periods  T  of  the  slowest  oscillation. 

2.  Determination  of  the  periods  is  simpler  for  the  two  extreme  cases 
represented  by  ?1  =  0  and  5^  =  1  .  Consider  the  first  case,  5^  =  0  . 
Since  J^O)  =  0  ,  from  equation  E23 

tan  K1/2  =  0  (E24) 

and  hence 

K1/2  =  7T  ,  k  =  tt2 

Substitute  this  into  equation  E8  to  obtain 

T  -  ~B=  (E25) 

M> 

a  well  known  result. 

E5 


3.  Taking  the  second  case  ^  =  l  ,  as  J  (2K1/2)  is  finite, 
equation  E23  yields  C 


J2  (2K1/2)  = 


(E26) 


the  root  of  the  smallest  value  being 


2K1/2  =  3. 


which  gives 


3.8317 


K  =  4.0530 


Substituting  this  in  equation  E8,  one  has 


T  =  3.121 


(E27) 


4.  The  determination  of  the  roots  of  equation  E23  for  other 
values  of  ^  must  be  carried  out  by  trials.  Some  other  determinations 
of  K  are  given  in  table  El,  along  with  the  corresponding  values 
of  N  in  the  expression 


T  =  Bvdr’'"^ 


(E28) 


Table  El 


Values  of 

and  N  for  Various  Values  of 

K1/2 

N 

0.0 

3-142 

2.000 

0.10 

3.137 

2.002 

0.20 

3-100 

2.026 

0.30 

3.007 

2.090 

0.33 

2.952. 

2.128 

0.40 

/ 

2.837 

2.214 

0.50 

2.642 

2.378 

0.60 

2.454 

2.560 

0.70 

2.285 

2.750 

o.8o 

2.141 

2.934 

0.90 

2.020 

3.110 

0.75 

2.210 

2.843 

1.00 

2.013 

3.121 

